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GP METHOD SEPARATES SOLUTION REGIONS AND 
. UNCOVERS UNEXPECTED ANOMALIES 
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INCREASE IN P 0 PRODUCES DRAMATIC WEIGHT 
INCREASE: BASELINE MISSION NOT REALIZED 
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OPTIMAL WALL THICKNESS VARIES LINEARLY 
WITH OPTIMAL BUMPER THICKNESS 
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OPTIMAL DESIGN IS SENSITIVE TO BUMPER/WALL SEPARATION 
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PROVIDES EFFECT OF THREAT ON OPTIMAL DESIGN 



PROJECTILE DIAMETER (CM) 
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OPTIMAL DESIGN IS SENSITIVE TO THREAT 
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GP METHOD PROVIDES THE MINIMUM BUMPER/WALL 
SEPARATION REQUIRED FOR VARIOUS THREAT SCENARIOS 

1 PROJECTILE 

„T 40 20 10 5 1 VELOCITY 


r 



PROJECTILE DIAMETER (d) (cm) 


u S ^ 
O V u 
>-3 ^ m 
M 2 

> Q 2 

w z ^ 
d *< 


O *2 H 

§ 5 §d 

to fc O 

^ 55 a 

H 05 > 
U M 
b H M 
b M H 
b 3 B 

SSS 

£ b « 
-*>-3 33 

| P p 

* w 5 

CO O • 

H 05 £ 
Mbs 
co O 
_ _ co 55 

0 O < 

g g 5 

H <3 b 
H > M 
M * M 
05 o J 
DC b ^ 

*- z £ 

* s i 

z a | 

M ” 2 
DC Z D 
H O CQ 


OPTIMAL DESIGN IS NOT SENSITIVE TO PROJECTILE 
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OPTIMAL DESIGN IS NOT SENSITIVE TO PROJECTILE 
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OPTIMAL DESIGN IS HIGHLY SENSITIVE TO MISSION RISK 
FOR Po ABOVE 0.97 
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OPTIMAL DESIGN IS HIGHLY SENSITIVE TO MISSION DURATION IN 
* 10-30 YEAR REGION 
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DESIGN ERROR INDUCED BY REDUCING PROBLEM 
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BASELINE DESIGN PARAMETERS 
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DEBRIS ENVIRONMENT DRIVES FOR 
BASELINE CORE MODULE AREAS 
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TOTAL CORE MODULE SYSTEM AREA (m 
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METHOD CONFIRMS THE MINIMUM WEIGHT 
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OPTIMAL DESIGN IS SENSITIVE TO BUMPER/WALL SEPARATION 



BUMPER/WALL SEPARATION (cm) 
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OPTIMAL DESIGN IS NOT SENSITIVE TO PROJECTILE 
VELOCITY IN HIGH VELOCITY REGIONS 



PROJECTILE VELOCITY (km/sec) 






OPTIMAL DESIGN IS NOT SENSITIVE TO PROJECTILE 
VELOCITY IN HIGH VELOCITY REGIONS 






I 


OPTIMAL DESIGN IS HIGHLY SENSITIVE TO MISSION RISK 
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OPTIMAL DESIGN IS HIGHLY SENSITIVE TO MISSION 
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OPTIMAL DESIGN IS SENSITIVE TO MISSION DURATION IN 

10-30 YEAR REGION 
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WHAT YOU WILL SEE IN SECTION IV 
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THE BURCH PREDICTOR (NORMAL IMPACT) 
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BASELINE DESIGN PARAMETERS 
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BASELINE DESIGN PARAMETERS 
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R MODIFIED BURCH 
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DEBRIS ENVIRONMENT DRIVES FOR WILKINSON 
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VELOCITY SENSITIVITY FOR THE BOEING SUBPREDICTORS 
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SEPARATION BETWEEN BUMPER AND WALL (cm) 
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SENSITIVITY OF OPTIMAL DESIGN TO THREAT 
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WILKINSON MORE SENSITIVE TO MISSION RISK 
THAN MODIFIED BURCH 
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DEFINES OPTIMAL DESIGN DISTRIBUTIONS 











WHAT YOU WILL SEE IN SECTION V 
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BASELINE DESIGN PARAMETERS 
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LARGE DESIGN PENALTY INDUCED BY INCREASE IN 
PROBABILITY OF NO PENETRATION: BOEING PREDICTOR 
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DESIGN IS HIGHLY SENSITIVE TO MISSION DURATION. 
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QP A PP L IED m NYSMITH 


Consider the Nysmith Equation: 
t, = 5.08 U9 ) , 

with inequality constraints 
t x ld<. 0.5 , ydS 1.0 

Note that if 
d>0.24/tV^ 2 , then 


fi 

d 


<S0.5=* 


% > 5.08V 278 (0.24/i V^ 2 ) 139 /[(0.5) 528 ^ 139 ] or , 
d 



Therefore, the third inequality constraint of Nysmith may be written: 
d £ 0.24/i V 0 ' 2 

Letting the idealized "weight" be denoted by w , one has: 
w ~ + h 

Substituting Nysmith, one has: 
w = + 5.08 V™d Z9i l(t™h 1 ' 39 ) , 


with inequality constraints: 

77 £ 1.0 , 5.08 V 27 ^ 1 ' 92 /^ 528 ^ 139 ) < 1.0 

a 






A. Unconstrained Optimization 
The dual GP problem may be written: 
maximize: 
subject to: 

where: 

This is easily solved, giving: 

8, =0.35 , 82 = 0.65, 


v(8) = (1/8,)* 1 (Cj/82)^ 

8, + 8j = 1 

8,-0.52882 = 0 
c, = 5.08 V m d 2n lh 139 


w 0 = 5.52V 1 W 91 


r, = 1 . 93 VV* 1 /* 0 - 91 


= 3.59V lt d l ' 9l /h 091 


is the minimum "weight". 

is the optimal bumper thickness, and 

is the optimal wall thickness. 


Note that the optimal ratios are given by: 


* 1# =0.35 = 8, , —^— = 0.65 = 82 

'lo + '2o 


Furthermore, this analytic solution is valid for 


d <0.23/1 v ° 2 why? 


h > A35dV => r, o < 1.93V 18 d L91 (4.35dV°T or 


t, £0.5d , as desired. Also, 

l 0 

^=1.86 t lo <d , as required. 




Thus the constrained problem need only be solved for 


0.23/t V 02 <>d< 0.24/t V" 02 

B. Constrained Optimization 

The constrained dual GP problem may be written: 

maximize: (S A , , 2 A, A/, A» 

v(8) = (l/8 1 ) (c/Sj) (c 3 ) (c 4 ) 

subject to: 8, + Sj = 1 

8, - 0.52882 + 8,,' - 0.5288, 2 ' = 0 
where: c, is defined as in part A 

c 2 = c,/d 
c 3 = 2/d 

This is a 2 degree-of-difficulty problem, since there are 2 equations and 4 unknowns. It is solved 
numerically. 





DERIVATION OF CORE MODULE CONFIGURATION WEIGHT 
Let: r 12 = radius of a core module 

r 22 = r i2 + t 2 

r n = r 22 +h 

r 2i = r n + r i 

L = core module length 

Then, 

w=nL [p,(4 - rf,) + p 2 (4 - rfj)] 
is the total core module weight. 




Let: r 12 = 210.8 cm 


L = 1356.7 cm 


Pi = P 2 = 2-81 gm/cm 3 


Applying the Nysmith equation, one has: 

w CMC = 12 tl + c x t x + 122 c/ x 412 + 25700 c 2 / 1 ”°' 528 + 310c 2 V" 


where: 


Ci = 5059 + 24 h , and 


v ^d zn /h 139 


The chosen approximate weight is given by: 


w cMc Afpnt — c ih + 25700c 2 /! 


The unconstrained GP problem is then: 


maximize: 


v(8) = (c 1 /8 1 ) 


8, f 25700 c 2 


subject to: 81 + 82=1 

8i-. 52882 = 0 


This is easily solved yielding: 
81 =0.346 , 82 = 0.654 










WILKINSON 


Recall that: 


■<£!=> t 2 = 0.364D 3 p V(^ 2 Pz) and. 


> 1 =» V = 0.364O 4 P ; V^Vi'iPz) 

Pin 


Consider the case where: 


Then, w = fj + c,^ -1 is the idealized weight. 


Where: 


Cj =0.364D 4 p 3 V jV /(L 2 5 2 p 1 p 2 ) 

The dual GP problem may be written: 


maximize: A , * 

v(5) = (1/8,) (Cj/Sj) 

subject to: 8 1 + 8 2 =1 

5 1 -5 2 = 0 


Thus, 8 t = 8j = 1/2 , and 





Thus, the procedure is the following: 


1) Determine r lo , from the above equation. 

D p_ 

2) Compute — - . 

Pl'l 0 


D p. 

3) If > 1 , then quit. Optimal design is 

Pih a 


Dp. 

4) If TT“ ^ I » optimal design is 


(*!<> ’ r2 o] ’ 


Pi'i # 


f 


'■.•V 

\ 

l plr, JJ 


Optimal Ratios for Wilkinson 


Recall that the optimal thickness distribution between bumper and wall is equal unless: 



. or 


r, ^ 

l 0 


D P, 
Pi 


This implies: 


0.604 





» or 


D S> 1 .66s 


^ ^2p2 

V„P' 


\l/2 

/ 





GP APPLIED TO MODIFIED BURCH 


Recall the Modified Burch Predictor 
.where 


rJ*"’' 

\N 


fr'SP 9 

— /S 0 ' 71 , and 

vi 


.-0.57 


r = 2.8(r 1 /D) +1.58(r 1 /D)' 

Thus, 

w=r,+ 2.$CD~ OS7 t° SI + is the idealized weight. 


Thus, the GP dual problem may be written: 
Maximize: / , \ 8 i 


v(8) = 


'l Yf2.8rzr" 




J 


1.58£ D 0 57 ^ 


Subject to: 8 t + 0.5752 - 0.5783 = 0 , and 

81 + 82+83 = 1 


Thus; we have: 

52 = 2.33(1-1.5783) , 8 , = 1.33(25j-l) , 

and since 8 j > 0 , 83 > 0 , we have 0.5 < 5 3 < 0.64 . 

This is a one degree-of-difficulty problem with the following procedure: 

1) Vary 83 from 0.5 to 0.64 to find the max v( 8 ) . 

2 ) Using the corresponding 83 , solve for 8 , , S 2 . 

3) Let r lo = 8 1 (max(v(5))) . 

4) Let £ 2o = max(v( 8 ))-r lo . 




OPTIMAL DESION ALGORITHM FOR BOEING PREDICTOR 


1. Compute optimal design for PEN4 Predictor, ^r 1# , j . 

2. Check against PEN4 constraint, V < V^r 1# ID j+ 4000 ? 

3. If satisfied, the optimal design is = ^r 1# , ^ j . 

4. Else, compute optimal designs for Modified Burch, fr, # , 
Wilkinson, ^ 


, and 


Predictors. 


5. Compute Wilkinson wall induced by optimal Modified Burch bumper, r 2 


6. Compute Burch wall induced by optimal Wilkinson bumper, t 2 t x 


Find , %,) - nnn[(< s , max^ , ^ J . , 


f 


max 


* 2 * ’ X \ 




'//J 










METEOR 1 


METEOR 1 is the meteoroid environment model (see Survivability and 
Predictor Reference 1) used in the design of protective systems for 
spacecraft. It is used to define projectile mass/diameter as a preprocessor 
to IMPACT5 and its derivatives. 

METEOR 1 takes as inputs the spacecraft exposed surface area, the mission 
duration, probability of no penetration and altitude, and the particle 
density. METEOR 1 accounts for Earth shielding and gravitational 
defocussing factors. The output is the critical projectile diameter. A 
sample input (MET. IN), output (MET.OUT) and program listing 

(METEOR 1. LIS) follow. 




•NO. OF CASES 

! SURFACE AREA IN SQ. METERS 
DURATION IN YEARS 
! PROB . OF NO PENETRATION 
I ALT . IN KM 

I DENSITY IN GM/CU8IC CM 



INPUT 


SURFACE AREA IN SQUARE METERS = 400.0000 

MISSION DURATION IN YEARS = 10.00000 

PROB. OF NO PENETRATION = 0.9700000 

ALTITUDE IN KILOMETERS = 500.0000 

PROJ. DENSITY IN GM/CUBIC CM = 0.5000000 

OUTPUT 


DIAM= 0.4614225 



12-Aug-1987 13:47:19 VAX FORTRAI^ 

12-Aug-1987 13:47:07 SAI USRDISK 



I 

I 

0 

I 

t 

0 

I 

0 

I 

0 

I 

0 

I 

01 

B: 

Ol 

I 

O' 

B 

01 

I 

J 


pool 
0002 
0003 
10004 
1)005 
0006 
|)007 
10008 
0009 
*0010 
loon 
1)012 
0013 
|)014 
P015 
0016 
J0017 
IDO 18 
*0019 
0020 
1)021 
1)022 
0023 
J0024 
1)025 
1)026 
0027 
1)028 
1)029 
0030 
W)031 
1)032 
1)033 
0034 
1)035 
1)036 
0037 
W)038 
1)039 
1)040 
0041 
1)042 
f0O43 
k0044 
p04b 


OPEN (UNIT=9, TYPE® 'OLD' , ACCESS 3 1 SEQUENTIAL 1 ,NAME =, MET.IN' ) 
0PEN(UNIT=10,TYPE='NEW' .ACCESS®' SEQUENTIAL' ,NAME='MET. OUT' ) 
READ(9,*)NCASES 
DO 10 1-1 ,NCASES 
READ(9,*)SA 
READ(9,*)T 
READ(9,*)P0 
READ(9,*)ALT 
READ(9,*)DENS 
WRITE (10,*)' INPUT' 

WRITE (10,*) 

WRITE(10,*)' SURFACE AREA IN SQUARE METERS = ' ,SA 

WRITE(10,*)' MISSION DURATION IN YEARS = ' ,T 

WRITE(10,*)' PROB. OF NO PENETRATION = ' ,P0 

WRITE(10,*)' ALTITUDE IN KILOMETERS = ' ,ALT 

WR I TE ( 10 , * ) ' PROJ. DENSITY IN GM/CUBIC CM = '.DENS 

WRITE (10,*) 

WRITE (10,*)' OUTPUT' 

WRITE (10,*) 

T=31536000.*T 

FLUX=-1 ,*ALOG(PO)/(SA*T) 

RA=6371 ./ (6371. +ALT) 

G£=.568+.432*RA 

THETA=ATAN(6371 ./SQRT(ALT*( ALT+2 .*6371 . ) ) ) 

S=( 1 ,+COS(THETA) )/2. 

FLUX=FLUX/(GE*S) 

F=AL0G10( FLUX) 

IF(F.GE.-4.403)THEN 
WRITE (10,*)' MASS IS TOO SMALL' 

GO TO 10 
ENDIF 

IF(F.GT.-7.103.AND.F.LT.-4.403)THEN 

RAD=2.509-.25*(14.339+L) 

XM=10.**( (-1 .584+SQRT(RAD) )/ .125) 

ENDIF 

I F ( F .LE .-7 .103. AND.F.GE.- 14.37 ) THEN 
XM=10.**( (14.37+F)/-1 .213) 

ENDIF 

IF(F.LT.-14.37)THEN 
WRITE ( 10 ,*) ' MASS IS TOO LARGE' 

GO TO 10 
ENDIF 

D=(1.91*XM/DENS)**.333 
WRITE( 10,*) ' DIAM=',D 

CONTINUE 
END 



OF. POOR QUALITX 



debrisi 


DEBRIS 1 is the space debris environment model (see Survivability and 
Predictor Reference 2) in the design of protective systems for spacecraft. 
It is used to define projectile mass/diameter as a preprocessor to IMPACT5 
and its derivatives. 

DEBRISI takes as inputs the spacecraft projected debris area, the mission 
duration, probability of no penetration, and altitude (currently fixed at 500 
Km). The output is the critical projectile diameter. A sample input 
(DEB.IN), output (DEB. OUT), and program listing (DEBRISI. LIS) follow. 




! NO. OF CASES 

! PROJECTED AREA IN SQ. METERS 
! DURATION IN YEARS 
1PROB. OF NO PENETRATION 



INPUT 


PROJ. AREA IN SQUARE METERS 
MISSION DURATION IN YEARS = 
PROB. OF NO PENETRATION = 

OUTPUT 


574.0000 

10.00000 

.9700000 


DIAM= 0.8446254 




I 


12-Aug-1987 13:52:31 
12-Aug-1987 13:52:24 


1 

0001 

0PEN(UNIT a 9 ,TYPE a ' OLD ' ,ACCESS a ' SEQUENTIAL ' ,NAME a ' 

1 

0002 

0PEN(UNIT a 10,TYPE a ' NEW' ,ACCESS=' SEQUENTIAL' ,NAME a 


0003 

READ(9,*)NCASES 

■ 

0004 

DO 10 I a l .NCASES 

1 

0005 

READ(9,*)AP 

0006 

READ(9,*)T 

_0007 

READ(9,*)P0 

WRITEdO,*)' INPUT 1 

■ 

0008 

■ 

0009 

WRITE(10 ,*) 


0010 

WRITEdO,*)' PROJ. AREA IN SQUARE METERS 

■ 

0011 

WRITEdO,*)' MISSION DURATION IN YEARS a 

| 

0012 

WRITEdO,*)' PROB. OF NO PENETRATION a 1 


0013 

WRITE(10,*) 


0014 

WRITEdO,*)' OUTPUT' 

1 

b015 

WRITEdO,*) 

■ 

0016 

FLUX=-1.*AL0G(P0)/(AP*T) 


0017 

F=AL0G10(FLUX) 

■0018 
po 19 

IF(F.GE.-5.46)THEN 

D=10.**( (F+5.46)/-2.52) 

0020 

ENDIF 

—0021 

IF(F.LE.-7.0)THEN 

B 

0022 

D=10.**( (F-21 .67 )/-10.32) 

■ 

0023 

ENDIF 

0024 

IF(F.LT.-5.46.AND.F.GT.-7.0)THEN 

■ 

0025 

D=10 .**( (5 .46+F)/- .63) 

10026 

ENDIF 

0027 

WRITE(10,*) ' DIAM= ' ,D 


JO 28 10 

CONTINUE 


J029 

END 


,P0 


VAX FORTRAN 
SAI USRDISK 
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OF POOR QUALITY 



IMPACT? 


IMPACT5 is a spacecraft protective systems design optimization code. 
IMPACT5 employs the Geometric Programming optimization technique in 
evaluating a number of nonlinear piecewise continuous, functional impact 
predictors. These include the Nysmith, Boeing, Madden, Wilkinson, Burch, 
and PEN4 predictors. Inputs vary depending on the predictor used, 
however, typical inputs include projectile characteristics (as determined 
from METEOR 1 or DEBRIS 1), design material properties, and general design 
configuration. In particular, IMPACT5 is an optimization code for a single 
bumper, single wall configuration. Outputs include the optimal design 
thicknesses (bumper and wall) and the minimum design weight. Sample 
input (IMPACT5.IN), output (IMPACT5.0UT) and program listing 
(IMPACT5.LIS) follow. 





6 

! NUMBER OF CASES 

1 

! NYSMITH PREDICTOR 

10. 

!PR0J. VELOCITY IN KM/SEC 

0.84 

!PR0J. DIAMETER IN CM 

10. 

! BUMPER/WALL SEPARATION IN CM 

2 

! BOEING PREDICTOR 

10. 

! PROJ . VELOCITY IN KM/SEC 

0.84 

! PRO J . DIA. IN CM 

2.81 

! PROJ . DENSITY IN GM/CUBIC CM 

2.81 

! BUMPER DENSITY IN GM/CUBIC CM 

2.81 

! WALL DENSITY IN GM/CUBIC CM 

10. 

! BUMPER/WALL SEPARATION IN CM 

0.401 

!WALL MATERIAL CONSTANT 

7344000. 

! BUMPER YIELD STRENGTH IN LB/SQUARE FOOT 

7344000. 

! WALL YIELD STRENGTH IN LB/SQUARE FOOT 

0 . 

! IMPACT ANGLE FROM NORMAL 

0.85 

! NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 

7.239E11 

! YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 

3 

! MADDEN PREDICTOR 

10. 

! PROJ . VELOCITY IN KM/SEC 

0.84 

1 PROJ . DIA. IN CM 

2.81 

! PROJ . DENSITY IN GM/CUBIC CM 

10. 

1BUMPER/WALL SEPARATION IN CM 

2.81 

! BUMPER/WALL DENSITY IN GM/CUBIC CM 

4 

1WILKINS0N PREDICTOR 

10. 

! PROJ • VELOCITY IN KM/SEC 

0.84 

! PROJ . DIA. IN CM 

2.81 

! PROJ . DENSITY IN GM/CUBIC CM 

2.81 

! BUMPER DENSITY IN GM/CUBIC CM 

2.81 

!WALL DENSITY IN GM/CUBIC CM 

10. 

! BUMPER/WALL SEPARATION IN CM 

0.401 

! WALL MATERIAL CONSTANT 

5 

! BURCH PREDICTOR 

10. 

I PROJ . VELOCITY IN KM/SEC 

0.84 

! PROJ . DIA. IN CM 

2.81 

! BUMPER DENSITY IN GM/CUBIC CM 

10. 

! BUMPER/ WALL SEPARATION IN CM 

0 . 

! IMPACT ANGLE FROM NORMAL 

0.85 

! NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 

7.239E11 

! YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 

6 

1PEN4 PREDICTOR 

1 . 

! PROJ . VELOCITY IN KM/SEC 

0.84 

IPROU. DIA. IN CM 

2.81 

! PROJ . DENSITY IN GM/CUBIC CM 

2.81 

! BUMPER DENSITY IN GM/CUBIC CM 

2.81 

!WALL DENSITY IN GM/CUBIC CM 

7344000. 

! BUMPER YIELD STRENGTH IN LB/SQUARE FOOT 

7344000. 

! WALL YIELD STRENGTH IN LB/SQUARE FOOT 

0 . 

! IMPACT ANGLE FROM NORMAL 




NYSMITH 


INPUT 

PROJECTILE VELOCITY IN KM/SEC = 10.00000 

PROJECTILE DIAMETER IN CM = 0.8400000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

OUTPUT 

BUMPER THICKNESS = 0.2575910 CM 

WALL THICKNESS = 0.4791193 CM 

MIN. WEIGHT = 0.7367103 CM 

CMC MIN. WEIGHT = 3799.036 KG 


BOEING 

INPUT 

PROJECTILE VELOCITY IN KM/SEC = 10.00000 

PROJECTILE DIAMETER IN CM * 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

BUMPER YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

WALL YIELD STRENGTH IN LB/SQUARE FT * 7344000. 

IMPACT- ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 
NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7.2389997E+11 

OUTPUT 

BUMPER THICKNESS = 9 .2028841E-02CM 

WALL THICKNESS = 0.4910776 CM 

MIN. WEIGHT = 0.5831065 CM 

CMC MIN. WEIGHT = 2978.630 KG 


MADDEN 

INPUT 

PROJECTILE VELOCITY IN KM/SEC = 10.00000 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM * 2.810000 

BUMPER/WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

OUTPUT 

BUMPER THICKNESS = 0.3788331 CM 

WALL THICKNESS = 0.3788331 CM 

MIN. WEIGHT = 0.7576661 CM 

CMC MIN. WEIGHT = 3934.167 KG 



WILKINSON 

INPUT 

PROJECTILE VELOCITY IN KM/SEC = 10.00000 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/ CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

OUTPUT 

BUMPER THICKNESS = 0.2128253 CM 

WALL THICKNESS = 0.2128253 CM 

MIN. WEIGHT = 0.4256506 CM 

CMC MIN. WEIGHT = 2208.558 KG 


MODIFIED BURCH 
INPUT 

PROJECTILE VELOCITY IN KM/SEC = 10.00000 

PROJECTILE DIAMETER IN CM = 0.8400000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

IMPACT ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 

NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7.2389997E+11 

OUTPUT 

BUMPER THICKNESS = 9 .2028841E-02CM 

WALL THICKNESS = 0.2539445 CM 

MIN. WEIGHT = 0.3459734 CM 

CMC MIN. WEIGHT = 1775.339 KG 


PEN4 

INPUT 

PROJECTILE VELOCITY IN KM/ SEC = 1.000000 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

WALL YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

IMPACT ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 

OUTPUT 

BUMPER THICKNESS = 2.1661645E-02CM 


WALL THICKNESS = 0 .OOOOOOOE+OOCM 

MIN. WEIGHT = 2.1661645E-02CM 

CMC MIN. WEIGHT = 222.9546 KG 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 

0012 25 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 35 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


0PEN( UNI T*10, TYPE* 'OLD' .NAME*' IMPACTS. IN' .ACCESS*' SEQUENTIAL' ) 
0PEN(UNIT*11, TYPE* 'NEW .NAME*' IMPACT5. OUT' .ACCESS* ' SEQUENTIAL ' ) 
READ(10,*)NRUNS 
DO 10 1*1 .NRUNS 
READ(10,*)NC0DE 
IF(NC0DE.EQ.1)G0 TO 25 
IF(NC0DE.EQ.2)G0 TO 35 
IF(NC0DE.EQ.3)G0 TO 45 
IF(NC0DE.EQ.4)G0 TO 55 
IF(NC0DE.EQ.5)G0 TO 65 
IF(NC0DE.EQ.6)G0 TO 75 
READ(10,*)V 
READ(10,*)D 
R£AD(10,*)H 

WRITE(11 ,*) ' NYSMITH' 

WRITE(11,*) 

WRITE (11,*) ' INPUT' 

WRITE(11,*) 

WRITE(11,*)' PROJECTILE VELOCITY IN KM/SEC = '.V 

WRITE(11 ,*) ' PROJECTILE DIAMETER IN CM = ' ,D 

WRITE(11,*)' BUMPER/WALL SEPARATION IN CM = ' ,H 

WR I TE (11 ,*) 

CALL N YSM I TH ( V , D , H ,T1 , T2 , WT , WTCMC ) 

WR I TE (11,*) ' OUTPUT’ 

WRITE(11,*) 

WRITE(11,*)' BUMPER THICKNESS = ' ,T1 , ' CM ' 

WRITE (11,*) ' WALL THICKNESS = * ,T2, ' CM' 

WRITE(11,*)' MIN. WEIGHT = '.WT.'CM' 

WRITE(11,*)' CMC MIN. WEIGHT * '.WTCMC, 'KG' 

WRITE(11,*) 

WRITE(11,*) 

WRITE(11,*) 

GO TO 10 

READ(10,*)V 

READ(10,*)D 

READ(10,*)RH0P 

READ(10,*)RH01 

READ(10,*)RH02 

READ(10,*)S 

READ(10,*)XL2 

READ(10,*)SY1 

READ( 10,*)SY2 

READ(10,*)THETA 

READ(10,*)XN 

READ(10,*)E1 

WRITE(11,*)' BOEING' 

WRITE(11,*) 

WRITE(1 1 ,*) ' INPUT' 

WRITE(11,*) 

WRITE(11,*)' PROJECTILE VELOCITY IN KM/SEC = ' ,V 

WRITE(ll t *)' PROJECTILE DIAMETER IN CM = ' ,D 

WRITE(11 ,*) ' PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 

WRITE(11,*)' BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 

WRITE(11,*)' WALL DENSITY IN GM/CUBIC CM = ' ,RH02 

WRITE(11,*)' BUMPER/WALL SEPARATION IN CM = ’ ,S 

WRITE(11,*)' WALL MATERIAL CONSTANT = ' ,XL2 

WRITE(11,*)' BUMPER YIELD STRENGTH IN LB/SQUARE FT * ' ,SY1 
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0058 


WRITEdl.*)' 

WALL YIELD STRENGTH IN LB/SQUARE FT * ' ,SY2 

0059 


WRITEdl.*) ' 

IMPACT ANGLE FROM NORMAL IN DEGREES * ' .THETA 

0060 


WRITEdl,*)' 

NUMBER OF PLATES TO PENETRATE AFTER FIRST', 

0061 

& 

' BUMPER = ’ ,XN 

0062 


WRITE(11,*)' 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE' , 

0063 

& 


' SEC = ' ,E1 

0064 


WRITEdl ,*) 


0065 


CALL B0EING(V,D 

,RH0P,RH01,RH02,S,XL2,SY1,SY2, THETA, 

0066 

& 

XN, 

El ,T1 ,T2 ,WT ,WTCMC ) 

0067 


WRITEdl,*) ' 

OUTPUT' 

0068 


WRITE(11,*) 


0069 


WRITEdl,*)' 

BUMPER THICKNESS * ' ,T1 , ' CM' 

0070 


WRITE(11,*)' 

WALL THICKNESS = ',T2,'CM' 

0071 


WRITE(11,*)' 

MIN. WEIGHT = ' ,WT ,' CM' 

0072 


WRITE(11,*)' 

CMC MIN. WEIGHT * ',WTCMC,'KG' 

0073 


WRITE(11,*) 


0074 


WRITE(11,*) 


0075 


WRITE(11,*) 


0076 


GO TO 10 


0077 

45 

REA0(10,*)V 


0078 


READ (10,* )D 


0079 


READ(10,*)RH0P 


0080 


READ dO,* )S 


0081 


REA0(10,*)RH0 

WRITE(11,*)' 


0082 


MADDEN' 

0083 


WRITE(11,*) 


0084 


WRITE (11,*) ' 

INPUT' 

0085 


WRITE(11,*) 


0086 


WRITEdl ,*) ' 

PROJECTILE VELOCITY IN KM/SEC = ' ,V 

0087 


WRITE(11,*)' 

PROJECTILE DIAMETER IN CM * ' ,D 

0088 


WRITE(11,*)' 

PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 

0089 


WRITE(11,*)’ 

BUMPER/WALL DENSITY IN GM/CUBIC CM = ’ ,RHO 
BUMPER/WALL SEPARATION IN CM * ' ,S 

0090 


WRITE(11,*)' 

0091 


WRITE(11,*) 

, R HOP , S , RHO , T1 , T2 , WT , WTCMC ) 

0092 


CALL MADDEN(V,D 

0093 


WRITE(11,*)' 

OUTPUT’ 

0094 


WRITE(11,*) 


0095 


WRITE(11,*) ' 

BUMPER THICKNESS = ' ,T1,'CM' 

0096 


WRITE{11,*)' 

WALL THICKNESS = ',T2,'CM' 

0097 


WRITE(11,*)' 

MIN. WEIGHT * ' ,WT,'CM' 

0098 


WRITE(11,*)' 

CMC MIN. WEIGHT = '.WTCMC, 'KG' 

0099 


WRITEdl,*) 


0100 


WRITE ( 11 ,* ) 


0101 


WRITE(11,*) 


0102 


GO TO 10 


0103 

55 

REA0{10,*)V 


0104 


READ(10,*)D 


0105 


READ(10,*)RH0P 

original; page is 

0106 


READ(10,*)RH01 

OF POOR QUALITY 

0107 


READ(10,*)RH02 

0108 


REA0(10,*)S 


0109 


READ(10,*)XL2 

WRITEdl,*)' 


0110 


WILKINSON' 

0111 


WRITE{11,*) 


0112 


WRITE(11,*)' 

INPUT' 

0113 


WRITE{11,*) 


0114 


WRITE(11,*)' 

PROJECTILE VELOCITY IN KM/SEC = ' ,V 
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011b 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 
015b 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 
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PROJECTILE DIAMETER IN CM = ' ,D 
PROJECTILE DENSITY IN GM/CUBIC CM = 1 ,RHOP 
BUMPER DENSITY IN GM/CUBIC CM = 1 ,RH01 
WALL DENSITY IN GM/CUBIC CM = ' ,RH02 
BUMPER/ WALL SEPARATION IN CM = ' ,S 
WALL MATERIAL CONSTANT = ‘ ,XL2 


WRITEUl,*)' 

WRITEU1,*)' 

WRITEdl,*)' 

WR I TE (11,*) 1 
WRITE (11 ,*) ' 

WR I TE (11,*) 1 
WRITE (11,*) 

CALL WILKINSON( V,D,RH0P,RH01 ,RH02,S,XL2, 
T1,T2,WT,WTCMC) 


65 


WRITEdl,*)' 
WRITE (11,*) 
WRITEdl ,*) 1 
WRITE(11 ,*) 1 
WRITEdl,*)' 
WRITE(11 ,*) 1 
WRITE(11,*) 
WRITE(11,*) 
WRITE(11,*) 

GO TO 10 

READ(10,*)V 

READ(10,*)D 

READ(10,*)RH01 

READ(10,*)S 

READ(10,*)THETA 


OUTPUT' 

BUMPER THICKNESS = ' ,T1,'CM' 
WALL THICKNESS = ' ,T2 , ' CM ' 

MIN. WEIGHT = ' ,WT, ' CM' 

CMC MIN. WEIGHT = ' ,WTCMC,'KG' 



READdO,*)XN 

READ(10,*)E1 


i kirk ** 

MODIFIED BURCH ***** 


WRITE(11,*)' 

MODIFIED BURCH' 


WRITE(11,*) 

WRITE(11,*)' 

INPUT' 


WRITE(11,*) 

WRITE(11,*)' 

PROJECTILE VELOCITY IN KM/SEC = ' ,V 


WRITE(11,*)’ 

PROJECTILE DIAMETER IN CM = ' ,D 


WRITE(11,*)' 

BUMPER DENSITY IN GM/CUBIC CM * ' ,RH01 


WRITE(11,*)' 

BUMPER/WALL SEPARATION IN CM * ' ,S 


WRITE(11,*)’ 

IMPACT ANGLE FROM NORMAL IN DEGREES = '.THETA 


WRITE(11,*)' 

NUMBER OF PLATES TO PENETRATE AFTER FIRST', 

& 


' BUMPER = ' ,XN 


WRITE(11,*)' 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE' , 

& 

' SEC 

- * ,E1 


WRITE(11,*) 

WRITE(11,*)' 

WRITE(11,*) 

OUTPUT' 


CALL BURCH(V,D,RH01,S, THETA, 

& 

XN,E1 

,T1 ,T2,WT,WTCMC) 


WRITE(11,*)’ 

BUMPER THICKNESS = ' ,T1 , ' CM' 


WRITE(11,*)’ 

WALL THICKNESS = ' ,T2 , ' CM' 


WRITE(11,*)' 

MIN. WEIGHT = ' , WT , ' CM ' 


WRITE(11,*)' 
WRITE(11,*) 
WRITE(11,*) 
WRITE(11,*) 
GO TO 10 

CMC MIN. WEIGHT = ' ,WTCMC,'KG' 

75 

READ(10,*)V 
READ( 10,*)D 



READ(10,*)RH0P 



READ(10,*)RH01 
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0172 


READ(10,*)RH02 


0173 


READ(10,*)SY1 


0174 


READ(10,*)SY2 


0175 


READ(10,*)THETA 


0176 


WRITE(11,*) 1 

PEN4' 

0177 


WRITE(11 ,*) 


0178 


WRITE ( 11 ,*) 1 

INPUT' 

0179 


WRITE (11,*) 


0180 


WRITE(11,*)' 

PROJECTILE VELOCITY IN KM/SEC = ' ,V 

0181 


WRITE( 11 ,*) 1 

PROJECTILE DIAMETER IN CM = ' ,D 

0182 


WRITE(11,*)' 

PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 

0183 


WRITE (11 ** ) ' 

BUMPER DENSITY IN GM/CUBIC CM = 1 ,RH01 

0184 


WRITE(11,*)‘ 

WALL DENSITY IN GM/CUBIC CM = ' ,RH02 

0185 


WRITEdl ,*) ' 

BUMPER YIELD STRENGTH IN LB/SQUARE FT = ' ,SY1 

0186 


WRITE(1 1 ,*) 1 

WALL YIELD STRENGTH IN LB/SQUARE FT = ' ,SY2 

0187 


WRITE(11,*)‘ 

IMPACT ANGLE FROM NORMAL IN DEGREES = ' , THETA 

0188 


WRITE(11,*) 


0189 


CALL PEN4(V,D,RH0P 

,RH01 ,RH02,SY1, SY2, THETA, 

0190 

& 

T1,T2,WT,WTCMC) 

0191 


WRITE(11,*)' 

OUTPUT' 

0192 


WRITE(11,*) 


0193 


WRITE(11,*) 1 

BUMPER THICKNESS = ' ,T1 , ' CM* 

0194 


WRITE(11,*)' 

WALL THICKNESS * ' ,T2,'CM' 

0195 


WRITE(11,*)' 

MIN. WEIGHT = ' ,WT, ' CM' 

0196 


WRITE(11,*)' 

CMC MIN. WEIGHT * ' ,WTCMC,'KG' 

0197 


WRITE(11,*) 


0198 


GO TO 10 


0199 

10 

CONTINUE 


0200 


STOP 


0201 


END 
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0001 

SUBROUTINE NYSMITH(V,D,H,T1,T2,WT,WTCMC) 

0002 

DMAX=0 ,24*H*V**-0 .2 

0003 

IF(D.GT.DMAX)THEN 

WRITE(11,*)'N0 SOLUTION— PROJ. DIA. TOO 

0004 

0005 

GO TO 16 

0006 

ENOIF 

0007 

Tl=l .93*V**0 .18*D**1 .91/H**0.91 

0008 

T2=1.86*T1 

0009 

WT=T1+T2 

0010 

WTCMC=Tl**2.+2 .*T1*(211 .+T2+H)+T2**2 .+422 

0011 

WTCMC=12.*WTCMC 

0012 16 

CONTINUE 

0013 

RETURN 

0014 

END 


PROGRAM SECTIONS 


Name 

Bytes 

Attributes 







0 $CODE 

214 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

1 SPDATA 

45 

PIC CON REL 

LCL 

SHR 

NOEXE 

RO 

NOWRT 

LONG 

2 $ LOCAL 

8 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 

Total Space Allocated 

267 









ENTRY POINTS 

Address Type Name 
0-00000000 NYSMITH 


VARIABLES 


Address Type 

Name 

Address 

Type 

Name 

Address Type 

Name 

AP-000000080 R*4 

D 


R*4 

DMAX 

AP-OOOOOOOCO R*4 

H 

AP -000000140 R*4 

T2 

AP -000000040 

R*4 

V 

AP-000000180 R*4 

WT 


LABELS 

Address Label 
0-000000D5 16 
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0001 

0002 

0003 

0004 
UU05 
0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


SUBROOTINE BOEING( V ,D,RHOP ,RH01 ,RH02 ,S,XL2 ,SY1 ,SY2, THETA, 
& XN,E1,T1 ,T2,WT,WTCMC) 

***** p£N4 ***** 

Tl=0.16 
V=V* 3280. 

0=0/30.48 
RP=D/2 .0 
RH0P=RH0P*1 .94 
RH01=RH01*1 .94 
RH02»RH02*1 .94 
NITSP=0 
NITSP=NITSP+1 
NP1=0 

TlP=Tl/30.48 

T2P=FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIP) 
WT*T1P+T2P 
IF( NITSP • EQ . 1 ) THEN 
T1P1«1.1*T1P 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 
WT1-T1P1+T2P1 
ENDIF 

IF(WT1.GT.WT)THEN 

T1P1=0.82*T1P1 

T2P1»FT2P(RH0P,V,RP,SY1, THETA, RH01 ,SY2,D,RH02, TIPI) 
WT1=T1P1+T2P1 

590 IF(WT1.GT.WT)THEN 

GO TO 601 
ELSE 

T1P-T1P1 

T2P-T2P1 

WT=WT1 

NP1=NP1+1 

IF(NP1 .EQ. 100) THEN 
WRITE(11,*) ' NO CONVERGENCE IN PEN4' 

GO TO 557 
ENDIF 

T1P1=0.9*T1P1 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 
WT1=T1P1+T2P1 
GO TO 590 
ENDIF 
ELSE 

579 T1P-T1P1 

T2P=T2P1 
WT=WT1 

T1P1=1.1*T1P1 

T2P1=FT2P(RH0P , V ,RP ,$Y1 .THETA, RH01 ,SY2 ,D ,RH02 , TIPI ) 
WT1=T1P1+T2P1 
IF(WT1.GT.WT)THEN 
GO TO 601 
ELSE 

NP1=NP1+1 

IF(NP1.EQ. 100) THEN 
WRlTE( 11 ,*) * NO CONVERGENCE IN PEN4 ' 

GO TO 557 
ENDIF 
GO TO 579 


BOEING 
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0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 


ENDIF 
END IF 

601 CONTINUE 
D=30.48*D 
RH0P=RH0P/1 .94 
RH01=RH01/1 .94 
RH02=RH02/1 .94 
T1P=30.48*T1P 
T2P=30.48*T2P 
IF(T1P/D.LE.0.4)VF=4100 
IF(T1P/D.GT.0.4)VF=4986*(T1P/D)**0.21 
VF=VF+4000. 

IF(V.LE.VF) THEN 

WRITE(11,*) 'INSIDE OF PEN4 LIMITS' 
T1=T1P 
T2=T2P 
GO TO 1102 
ENDIF 

557 CONTINUE 


0077 ***** WILKINSON ***** 


0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 

0113 

0114 


V-V/ 3280. 

T1=0.604*D**2.*RH0P/S 
T1=T1*SQRT(V/(XL2*RH01*RH02) ) 

RATI0=D*RH0P/(T1*RH01) 

IF( RATIO. GT.1.0)T2*T1 
IF(RATI0.LE.1.0)T2=T1/RATI0 
***** MODIFIED BURCH ***** 

V B=V*3280. 

* DB a D/2.54 
CM=SQRT(E1/RH01) 

CM*CM/30.48 

SB=S/2.54 

IF(THETA.LE.0.001)G0 TO 125 
CHI=TAN( THETA) -0.5 

F2=0.5-1.87*(T1B/D)+(5.*T1B/D-1.6)*CHI**3.0 
F2=F2+(1.7-12.*T1B/D)*CHI 
F3=0.32*(T1B/D)**0 .83 

F3=F3+0.48*(T1B/D)**0.33*(SIN(THETA) )**3 .0 
T2F=D*( (F1+0.63*F2)/XN)*(CM/V)**2.29 
T2F=T2F*(D/S)**0 .71 
T2N=F3*(CM/V)**1.33*D/XN 
IF(T2N.GE.T2F)T2B=T2N 
IF(T2N.LT.T2F)T2B=T2F 
T2B=T2B*2 .54 
IF(T2B.GT.T2)NREGI0N=3 
IF(T2B.GT.T2)T2=T2B 
GO TO 155 
125 CONTINUE 

TlBl=Tl/2.54 

NITSB=0 

XK1=(DB/XN)**1 . 7 1*( CM/VB)**2 .29/ SB**0 .71 

VDELTA=0.0 

DELTA30.52 

1099 DELTA2*2 .33*( 1 .-1 .57*DELTA3) 

DELTA1=1.33*(2.*DELTA3-1.) 

VDELTAl® ( 1 ./DELTA1 ) **DELTA1* ( 2 .8*XK1/ (DELTA2*DB**0 . 57 ) ) **DELTA2 
VDELTA1=VDELTA1*(1.58*XK1*DB**0.57/DELTA3)**DELTA3 



BOEING 
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0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 
012b 
0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 
0143. 
0144 
014b 

0146 

0147 

0148 

0149 
OlbO 
0151 
01b2 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 
0163 


IF(VDELTA1.LT.VDELTA)THEN 
DELTA1*1 ,33*(2.*DELTA3-1 .04) 

T1B 3 DELTA1*VDELTA 
T2B 3 VDELTA-T1B 
GO TO 499 
ENDIF 

VDELTA*VDELTA1 
DELTA3 3 DELTA3+0 .02 
IF(DELTA3.GT.0.63)THEN 
T1B 3 DELTA1*VDELTA 
T2B-V0ELTA-T1B 
GO TO 499 
ENDIF 

GO TO 1099 
499 CONTINUE 

***** COMPARISON OF MODIFIED BURCH AND WILKINSON ***** 

199 CONTINUE 

T10W 3 Tl/2.54 

F10W 3 1.58*(DB/T10W)**0.57+2.80*(T10W/DB)**0.57 
T2BT10W 3 ( F10W/XN)**1 .71*(CM/VB)**2 .29*DB**1 .71 
T2BT10W 3 T2BT10W/SB**0 .71 
T2BT10W 3 T2BT10W*2 .54 
RATIOB 3 (DB*RHOP ) / ( RH01*T1B ) 

T2WT10B*0 .364*0**3 .*RHOP*V/ ( XL2*RH02*S**2 . ) 
IF(RATI0B.GT.1.0)T2WT10B 3 T2WT10B*RATI0B 
I F( T2BT10W.GT .T2 ) T2=T2BT10W 
T2B»T2B*2.54 

I F ( T2WT10B . GT.T2B ) T2B*T2WT10B 
T1B»T1B*2.54 • 

IF(T1B+T2B.LT.T1+T2)THEN 

T1 3 T1B 

T2=T2B 

ENDIF 

155 CONTINUE 

1102 IF(T2.LE.0.01)THEN 

WRITE(11,*) 'PRESSURE CONSTRAINT ON WALL IN EFFECT' 
SIGMA 3 SY2/144. 

T2=3099. 1/SIGMA 
ENDIF 
WT 3 T1+T2 
R12 3 211 . 

R22*211.+T2 
Rll 3 211.+T2+S 
R21*211 .+T1+T2+S 
VB 3 4.27*(R21**2.-R11**2.) 

VW 3 4.27*(R22**2.-R12**2.) 

WTCMC=RH01*VB+RH02*VW 

156 RETURN 
END 
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0001 FUNCTION FT2B(DB,T1B,XN,CM,VB,SB) 

0002 F1*2.42*(DB/T1B)**0 .33+4.26*(TlB/DB)**0.33 

0003 F1»F1-4.18 

0004 FT2B=(F1/XN)**1 .71*(CM/VB)**2.29*DB**1 .71/SB**0 .71 

0005 RETURN 

0006 END 


PROGRAM SECTIONS 


Name 

Bytes 

0 $CODE 

163 

2 $LOCAL 

4 

Total Space Allocated 

167 


Attributes 

PIC CON REL LCL SHR EXE RD NOWRT LONG 

PIC CON REL LCL NOSHR NOEXE RO WRT LONG 


ENTRY POINTS 


Address Type Name 
0-00000000 R*4 FT2B 

VARIABLES 


Address Type 

Name 

Address Type 

Name 

Address 

Type 

Name 

AP-000000100 R*4 

CM 

AP -000000040 R*4 

DB 

** 

R*4 

FI 

AP-0000U0080 R*4 

TIB 

AP-000000140 R*4 

VB 

AP-OOOOOOOCO 

R*4 

XN 
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0001 FUNCTION FT2P(RH0P,V,RP,SY1 , THETA, RH01,SY2,D,RH02, TIP) * 

0002 A=1.33*RH0P*(V*RP)**2. - 

0003 8=8.U*SY1*EXP(-3.125E-04*V) /COS (THETA) I 

0004 C=1 ,33*RH0P*RP**2 .0 ■ 

0005 D1=RP*RH01/C0S(THETA) 

0006 XK1»1.67*(RH0P/(2.*SY2))**0.31 ■ 

0007 XK1=XK1*(0.281*D*RH0P/RH02)**0.33 | 

0008 XK1=XK1*C0S( THETA) 

0009 C1P1=(A-B*T1P)/(C+D1*T1P ) _ 

0010 I F( Cl PI. LE. 0.001) THEN ■ 

0011 FT2P=0.0 ■ 

0012 GO TO 999 

0013 ENDIF ■ 

0014 FT2P=XK1*C1P1**0 .31 | 

0015 999 RETURN 

0016 END _ 


PROGRAM SECTIONS 


Name 

Bytes 

0 $CODE 

222 

2 SLOCAL 

4 

Total Space Allocated 

226 


Attributes 

PIC CON REL LCL SHR EXE RD NOWRT LONG 

PIC CON REL LCL NOSHR NOEXE RD WRT LONG 


ENTRY POINTS 


Address Type 
0-00000000 R*4 

Name 

FT2P 







VARIABLES 

Address Type 

Name 

Address Type 

Name 

Address Type 

Name 

** R*4 

A 

** 

R*4 

B 


R*4 

C 

AP-000000200 R*4 

D 

** 

R*4 

D1 

AP-000000180 

R*4 

RH01 

AP-00000004@ R*4 

RHOP 

AP-OOOOOOOC0 

R*4 

RP 

AP-000000100 

R*4 

SY1 

AP-000000280 R*4 

TIP 

AP-000000140 

R*4 

THETA 

AP-00000008@ 

R*4 

V 


LABELS 

Address Label 


U-000000D9 


999 
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U001 SUBROUTINE MADDEN(V,D,RH0P,S,RH0,T1 ,T2,WT,WTCMC) 

0002 V=V*1 .E05 

0003 T1=0 .009*SQRT ( V ) *RH0P*D**2 .0 

0004 T1=T1/ (S*RH0**1 .5) 

0005 T2=T1 

0006 WT=T1+T2 

0007 R12=211 . 

0008 R22=211.+T2 

0009 Rll=211 .+T2+S 

0010 R21=211 .+T1+T2+S 

0011 VB=4 .27*( R21**2 .-Rl 1**2 . ) 

0012 VW=4 . 27* ( R22**2 . -Rl 2**2 . ) 

0013 WTCMC=RHO* ( VB+VW ) 

0014 RETURN 

0015 END 


PROGRAM SECTIONS 


Name 

Bytes 

Attributes 

0 $CODE 

159 

PIC CON REL LCL 

Total Space Allocated 

159 



ENTRY POINTS 


Address Type Name 
0-00000000 MADDEN 

VARIABLES 


Address 

Type 

Name 

Address Type 

Name 

<S> 

00 

o 

o 

o 

o 

o 

o 

o 

1 

Q_ 

< 

R*4 

D 


R*4 

Rll 


R*4 

R22 

AP-000000140 

R*4 

RH0 

AP-000000180 

R*4 

T1 

AP-0000001C@ 

R*4 

T2 

** 

R*4 

VW 

AP-000000200 

R*4 

WT 


FUNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTHSSQRT 


SHR EXE RD NOWRT LONG 


Address Type Name 

** R*4 R12 

AP-OOOOOOOCO R*4 RHOP 
AP-000000040 R*4 V 
AP-000000240 R*4 WTCMC 
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0001 SUBROUTINE WILKINS0N(V,D,RH0P,RHQ1,RH02,S,XL2, 

0002 & T1 ,T2,WT,WTCMC) 

0003 T1=0. 604*0**2. *RH0P/S 

0004 T1®T1*SQRT(V/(XL2*RH01*RH02) ) 

0005 RATI0=D*RH0P/(T1*RH01) 

0006 IF(RATI0.GT.1.0)T2=T1 

0007 IF( RATIO. L£.1.0)T2 S T1/RATI0 

0008 WT=T1+T2 

0009 R12=211 . 

0010 R22=211.+T2 

0011 Rll*211.+T2+S 

0012 R21*211.+T1+T2+S 

0013 VB*4. 27*(R2 1**2. -R1 1**2. ) 

0014 VW=4 . 2 7* ( R2 2**2 . -R 1 2**2 . ) 

0015 WTCMC=RH01*V8+RH02*VW 

0016 RETURN 

0017 END 




VAX FO 



PROGRAM SECTIONS 

Name Bytes Attributes 

0 $C0DE 181 PIC CON REL LCL SHR EXE RO NOWRT LONG 

Total Space Allocated 181 

ENTRY POINTS 

» 

Address Type Name 

0-00000000 WILKINSON I 


VARIABLES 


Address Type 

Name 

Address Type 

Name 

Address Type 

Name 

AP -000000080 

R*4 

D 

** R*4 

Rll 

** R*4 

R12 

** 

R*4 

R22 

** R*4 

RATIO 

AP-000000100 R*4 

RH01 

AP-OOOOOOOC0 

R*4 

RHOP 

AP -000000 180 R*4 

S 

AP-000000200 R*4 

T1 

AP -000000040 

R*4 

V 

** R*4 

VB 

** R*4 

VW 

AP-OOOOOO2C0 

R*4 

WTCMC 

AP -000000 1C0 R*4 

XL2 




FUNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTHSSQRT 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
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***** MODIFIED BURCH ***** 

SUBROUTINE BURCH(V,D,RH01,S, THETA, 

& XN,E1 ,T1 ,T2,WT,WTCMC) 

VB=V*3280. 

DB=D/2.54 

CM=SQRT(E1/RH01) 

CM=CM/30.48 

SB*S/2.54 

I F (THETA. LE. 0.001) GO TO 425 
CHI=TAN( THETA) -0.5 

F2 a 0.5-1.87*(T1B/D)+(5.*T1B/D-1.6)*CHI**3.0 
F2=F2+(1.7-12.*T1B/D)*CHI 
F3 a 0.32*(TlB/D)**0 .83 

F3 a F3+0.48*(TlB/D) **0.33* ( SI N( THETA) )**3 .0 
T2F=D*((F1+0.63*F2)/XN)*(CM/V)**2.29 
T2F a T2F*(D/S)**0.71 
T2N a F3*(CM/V)**l .33*D/XN 
IF(T2N.GE.T2F)T2B=T2N 
IF (T2N.LT. T2F)T2B=T2F 
T2B a T2B*2 .54 
IF(T2B.GT.T2)NREGI0N a 3 
IF(T2B.GT.T2)T2 a T2B 
GO TO 499 
425 CONTINUE 

NITSB=0 

XK1= ( DB/XN ) **1 . 7 1* ( CM/VB ) **2 .29/ SB**0 

VDELTA a 0.0 

DELTA3=0.52 

1099 DELTA2 a 2.33*(l .-1.57*DELTA3) 

DELTAl a 1.33*(2.*DELTA3-l.) 

VDELTA1*( 1 ./DELTA1)**DELTA1*(2.8*XK1/ (DELTA2*DB**0.57 ) )**DELTA2 
VDELTA1*VDELTA1*( 1 .58*XK1*DB**0 .57/DELTA3 ) **DELTA3 
IF(VDELTA1.LT.VDELTA)THEN 
DELTAl a l .33*(2.*DELTA3-1.04) 

T1*DELTA1*VDELTA 
T2=VDELTA-T1 
GO TO 499 
ENDIF 

VDELTA*VDELTA1 
DELTA3=DELTA3+0 .02 
IF(DELTA3.GT.0.63)THEN 
T1 3 DELTA1*VDELTA 
T2 a VDELTA-Tl 
GO TO 499 
ENDIF 

GO TO 1099 
499 CONTINUE 

T1=T1*2.54 
T2=T2*2.54 
WT=T1+T2 
R12=211 . 

R22=211 .+T2 
Rll=211 .+T2+S 
R21=211 .+T1+T2+S 
VB=4.27*(R21**2.-R11**2.) 

VW=4.27*(R22**2 .-Rl 2**2 • ) 

WTCMC=RH01*VB+2.81*VW 



VAX 

SAI 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 588 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 577 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


SUBROUTINE PEN4( V ,D ,RHOP ,RH01 ,RH02 ,SY1 ,SY2 , THETA, 
T1,T2,WT,WTCMC) 

Tl=0.16 

V*V*3280. 

D=D/30.48 
RP=D/2.0 
RH0P=RH0P*1 .94 
RH01=RH01*1 .94 
RH02=RH02*1 .94 
NITSP=0 
NITSP=NITSP+1 
NP1=0 

TlP=Tl/30.48 

T2P=FT2P ( RHOP , V , RP , SY 1 , THETA , RH01 , SY2 ,D , RH02 , T1 P ) 
WT=T1P+T2P 
I F( N I TSP • EQ . 1 ) THEN 
T1P1*1 .1*T1P 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01 ,SY2,D,RH02, TIPI) 
WT1*T1P1+T2P1 
ENDIF 

IF(WT1.GT.WT)THEN 

T1P1=0.82*T1P1 

T2P1=FT2P(RH0P , V ,RP ,SY1 , THETA, RH01 ,SY2 ,D ,RH02 ,T1P1 ) 
WT1=T1P1+T2P1 
IF(WT1.GT.WT)THEN 
GO TO 599 
ELSE 

T1P=T1P1 

T2P=T2P1 

WT=WT1 

NP1=NP1+1 

I F ( NP1 .EQ. 100) THEN 
WRITE(11 ,*) 'NO CONVERGENCE IN PEN4* 

GO TO 555 
ENDIF 

T1P1=0.9*T1P1 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01 ,SY2,D ,RH02, TIPI) 
WT1=T1P1+T2P1 
GO TO 588 
ENDIF 
ELSE 

T1P=T1P1 

T2P=T2P1 

WT=WT1 

T1P1=1.1*T1P1 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01 ,SY2,D,RH02, TIPI) 
WT1=T1P1+T2P1 
IF(WT1.GT.WT)THEN 
GO TO 599 
ELSE 

NP1=NP1+1 

I F ( NP 1 .EQ. 100) THEN 
WRITE( 11 ,*) ' NO CONVERGENCE IN PEN4' 

GO TO 555 
ENDIF 
GO TO 577 
ENDIF 


PEN4 


10-Aug-1987 11:29:15 VAX FOf 

10-Aug-1987 11:29:03 SAI_USI 


0058 


END IF 

0059 

599 

CONTINUE 

0060 


0*30.48*0 

0061 


T1P*30.48*T1P 

0062 


T2P*30.48*T2P 

0063 


IF(T1P/D.LE.0.4)VF*4100 

0064 


IF(T1P/D.GT.0.4)VF=4986*( 

0065 


VF=VF+4000. 

0066 


IF(V.GT.VF)THEN 

0067 


WRITE(11,*) 'OUTSIDE OF 

0068 


GO TO 1100 

0069 


END IF 

0070 


T1*T1P 

0071 


T2*T2P 

0072 

555 

CONTINUE 

0073 


WT*T1+T2 

0074 


R12*211. 

0075 


R22-211.+T2 

0076 


R11-211.+T2+10. 

0077 


R21*211 .+T1+T2+10 . 

0078 


VB»4.27*(R21**2.-R11**2.) 

0079 


VW*4 .27*(R22**2 .-R12**2 . ) 

0080 


WTCMC»RH01*VB+RH02*VW 

0081 

1100 

RETURN 

0082 


END 


PROGRAM SECTIONS 


Name 


Bytes 

Attributes 







0 SCODE 


775 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

1 SPDATA 


44 

PIC CON REL 

LCL 

SHR 

NOEXE 

RD 

NOWRT 

LONG 

2 $LOCAL 


120 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 

Total Space Allocated 


939 








ENTRY POINTS 










Address Type Name 










U-OOUOOUUO PEN4 











VARIABLES 

Address Type 

Name 

Address Type 

Name 

Address 

Type 

Name 

AP -000000080 

R*4 

D 


1*4 

NITSP 

** 

1*4 

NP1 

** 

R*4 

R12 

★* 

R*4 

R21 

irk 

R*4 

R22 

AP-000000140 

R*4 

RH02 

AP-OOOOOOOCO 

R*4 

RHOP 

2-00000000 

R*4 

RP 

AP-OOOOOOICO 

R*4 

SY2 

AP-000000240 

R*4 

T1 

2-00000004 

R*4 

TIP 

AP-000000280 

R*4 

T2 

** 

R*4 

T2P 

** 

R*4 

T2P1 

AP -000000040 

R*4 

V 

** 

R*4 

VB 

2-OOOOOOOC 

R*4 

VF 



IMPACTgV 


IMPACT5V is a spacecraft protective systems design optimization code 
similar to IMPACT5. IMPACT5V differs from IMPACT5 in that the optimal 
design is weighted according to the chosen space debris velocity 
probability distribution. IMPACT5V employs the Geometric Programming 
optimization technique in evaluating a number of nonlinear piecewise 
continuous, functional impact predictors. These include the Nysmith, 
Boeing, Madden, Wilkinson, and Burch predictors. Inputs vary depending 
on the predictor used, however, typical inputs include the space debris 
velocity distribution file, projectile characteristics (as determined from 
METEOR 1 or DEBRIS 1), design material properties, and general design 
configuration. In particular, IMPACT5V is an optimization code for a 
single bumper, single wall configuration. Outputs include the optimal 
design thicknesses (bumper and wall) and the minimum design weight. 
The velocity distribution files for 500 Km altitude and 30 degree 
inclination (500KM30DEG.DAT) and 60 degree inclination (500KM60 
DEG.DAT) follow. The file being used must be assigned to FOR012 before 
running IMPACT5V. Sample input (IMPACT5V.IN), output 
(IMPACT5V.OUT), and program listing (IMPACT5V.LIS) follow the velocity 
distribution files. 




1 

0.005 

2 

0.005 

3 

0.02 

4 

0.05 

5 

0.10 

6 

0.05 

7 

0.075 

8 

0.08 

9 

0.07 

10 

0.23 

11 

0.155 

12 

0.08 

13 

0.05 

14 

0.02 

18 

0.005 

16 

0.005 




1 

2 

3 

4 
b 
6 

7 

8 
y 

10 

11 

12 

13 

14 

15 

16 


0.01 

0.01 

0.01 

0.01 

0.03 

0.03 

0.01 

0.01 

0.005 

0.09 

0.09 

0.07 

0.28 

0.28 

0.065 

0.005 




5 

1 

0.84 

10 . 

2 

0.84 

2.81 

2.81 

2.81 

10 . 

0.401 

7344000. 

7344000. 

0 . 

0.85 

7.239E11 

3 

0.84 

2.81 

10 . 

2.81 

4 

0.84 

2.81 

2.81 

2.81 

10 . 

0.401 

5 

0.84 

2.81 

io.- 

0 . 

0.85 

7.239E11 


! NUMBER OF CASES 

1 NYSMITH PREDICTOR 

1PR0J. DIAMETER IN CM 

! BUMPER/WALL SEPARATION IN CM 

! BOEING PREDICTOR 

i PROJ . DIA. IN CM 

! PRO J . DENSITY IN GM/CUBIC CM 

! BUMPER DENSITY IN GM/CUBIC CM 

! WALL DENSITY IN GM/CUBIC CM 

1BUMPER/WALL SEPARATION IN CM 

!WALL MATERIAL CONSTANT 

! BUMPER YIELD STRENGTH IN LB/ SQUARE FOOT 

1WALL YIELD STRENGTH IN LB/SQUARE FOOT 

! IMPACT ANGLE FROM NORMAL 

! NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 

! YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 

! MADDEN PREDICTOR 

i PROJ . DIA. IN CM 

! PROJ . DENSITY IN GM/CUBIC CM 

! BUMPER/WALL SEPARATION IN CM 

! BUMPER/ WALL DENSITY IN GM/CUBIC CM 

! WILKINSON PREDICTOR 

I PROJ . DIA. IN CM 

1PR0J. DENSITY IN GM/CUBIC CM 

! BUMPER DENSITY IN GM/CUBIC CM 

1WALL DENSITY IN GM/CUBIC CM 

! BUMPER/WALL SEPARATION IN CM 

! WALL MATERIAL CONSTANT 

! BURCH PREDICTOR 

! PROJ . DIA. IN CM 

! BUMPER DENSITY IN GM/CUBIC CM 

! BUMPER/ WALL SEPARATION IN CM 

! IMPACT ANGLE FROM NORMAL 

! NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 
! YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 


ORIGINAL PAGE IS 
OF POOR QUALITY 




NY SMITH 


INPUT 

PROJECTILE DIAMETER IN CM * 0.8400000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

OUTPUT 

BUMPER THICKNESS * 0.2498978 CM 

WALL THICKNESS = 0.4648098 CM 

MIN. WEIGHT * 0.7147076 CM 

CMC MIN. WEIGHT * 3685.385 KG 


BOEING 

INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

BUMPER YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

WALL YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

IMPACT ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 
NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7 .2389997E+U 

INSIDE OF PEN4 LIMITS 
INSIDE OF PEN4 LIMITS 
OUTPUT 

BUMPER THICKNESS = 0.1605087 CM 

WALL THICKNESS = 0.5884032 CM 

MIN. WEIGHT = 0.7489119 CM 

CMC MIN. WEIGHT = 3837.381 KG 


MADDEN 

INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM * 2.810000 

BUMPER/WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

OUTPUT 

BUMPER THICKNESS * 0.3522651 CM 

WALL THICKNESS = 0.3522651 CM 

MIN. WEIGHT = 0.7045302 CM 

CMC MIN. WEIGHT * 3657.869 KG 



WILKINSON 


INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

OUTPUT 

BUMPER THICKNESS = 0.1978996 CM 

WALL THICKNESS = 0.1978996 CM 

MIN. WEIGHT = 0.3957993 CM 

CMC MIN. WEIGHT = 2053.512 KG 


MODIFIED BURCH 
INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

IMPACT ANGLE FROM NORMAL IN DEGREES = 0 .OOOOOOOE+OO 

NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7 .2389997E+11 

OUTPUT 

BUMPER THICKNESS = 0.1434212 CM 

WALL THICKNESS = 0.6418241 CM 

MIN. WEIGHT = 0.7852453 CM 

CMC MIN. WEIGHT = 4017.877 KG 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 24 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 25 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 26 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 35 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


DIMENSION XPV( 16) 

OPEN (UN IT* 12, TYPE* 'OLD' ,ACCESS=' SEQUENTIAL ' ) 

0PEN(UNIT=10 ,TYPE= ' OLD' .NAME*' IMPACT5V.IN' .ACCESS* ' SEQUENTIAL ' ) 
0PEN(UNIT*11 .TYPE* ' NEW' .NAME* ' IMPACT5V.0UT' ,ACCESS= 1 SEQUENTIAL ' ) 
DO 24 1*1,16 

READ (12,*) IV,XPV(IV) 

CONTINUE 
READ(10,*)NRUNS 
DO 10 I-l.NRUNS 
READ(10,*)NC0DE 
IF(NC0DE.EQ.1)G0 TO 25 
IF(NC0DE.EQ.2)G0 TO 35 
IF(NCODE.EQ.3)GO TO 45 
IF(NC0DE.EQ.4)G0 TO 55 
IF(NC0DE.EQ.5)G0 TO 65 
READ(10,*)D 
READ(10,*)H 

WRITE(11 ,*) ' NYSMITH' 

WRITE (11 ,*) 

WRITE (11 ,*) 1 INPUT' 

WRITE (11 ,*) 

WRITE(11 ,*) ' PROJECTILE DIAMETER IN CM = ' ,D 

WRITE(11,*)' BUMPER/WALL SEPARATION IN CM = ' ,H 

WR ITE (11,*) 

T1T=0.0 
T2T=0 .0 
DO 26 J=1 ,16 
V=FLOAT(J ) 

CALL NYSMITH( V,D,H,T1,T2,WT,WTCMC) 

T1T=T1T+T1*XPV( J ) 

T2T=T2T+T2*XPV( J) 

CONTINUE 

T1=T1T 

T2=T2T 

WT=T1+T2 

WTCMC=Tl**2.+2 .*T1*( 211 .+T2+H)+T2**2 .+422 .*T2 

WTCMC=12.*WTCMC 

WRITE(11 ,*) ' OUTPUT' 

WRITE(11 ,*) 

WRITE( 11 ,*) ' BUMPER THICKNESS = ',T1,'CM' 

WRITE( 11 ,*) ' WALL THICKNESS = ' ,T2, 'CM' 

WR I TE (11,*) ' MIN. WEIGHT = ',WT,'CM' 

WR ITE (11,*) ' CMC MIN. WEIGHT = \WTCMC,'KG' 

WRITE(11 ,*) 

WRITE (11 ,*) 

WRITE(11,*) 

GO TO 10 

READ(10,*)D 

READ(10,*)RH0P 

READ(10,*)RH01 

READ(10,*)RH02 

READ(10 ,*)S 

READ(10 ,*)XL2 

READ ( 10 ,*)SY1 

READ(10 ,*)SY2 

READ(10,*)THETA 

READ(10,*)XN 
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VAX 

SAI 


I 


FORT 
USRDISI 


58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

173 

174 

175 

176 

177 
)78 
)79 
>80 
>81 
>82 
>83 
>84 
)85 
>86 
)87 
>88 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 


36 


45 


& 

& 


& 


READ (10»*)E1 

WRITE( 11 ,*) 1 

WRITE (11 ,* ) 

WRITE (11,*) 1 

WRITE( 11 ,*) 

WRITE(11 ,*) 1 

WRITEdl,*)' 

WRITE(11,*)' 

WRITEdl,*)' 

WRITEdl,*)' 

WRITEdl,*)' 

WRITEdl,*)' 

WRITE(11,*)' 

WRITE(11,*)' 

WRITE(11,*)' 

WRITE(11,*)' 


BOEING' 

INPUT' 

PROJECTILE DIAMETER IN CM * ' ,D 
PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 
BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 
WALL DENSITY IN GM/CUBIC CM = ' ,RH02 
BUMPER/WALL SEPARATION IN CM = ' ,S 
WALL MATERIAL CONSTANT = ’ ,XL2 
BUMPER YIELD STRENGTH IN LB/SQUARE FT = ' ,SY1 
WALL YIELD STRENGTH IN LB/SQUARE FT = ' ,SY2 
IMPACT ANGLE FROM NORMAL IN DEGREES = ' , THETA 
NUMBER OF PLATES TO PENETRATE AFTER FIRST' , 

1 BUMPER = ' ,XN 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE' , 

SEC = ' ,E1 


WRITE(11,*) 

T1T S 0.0 
T2T=0.0 
DO 36 J-1,16 
V=FLOAT( J ) 

CALL BOE ING( V ,D ,RHOP ,RH01 ,RH02 ,S ,XL2 ,SY1 , SY2 , THETA , 
XN,E1,T1,T2,WT,WTCMC) 

T1T=T1T+XPV(J)*T1 

T2T=T2T+XPV ( J )*T2 

CONTINUE 

T1-T1T 

T2=T2T 

WT=T1+T2 

R12=211. 

R22-211.+T2 

Rll a 211.+T2+S 

R21 a 211 .+T1+T2+S 

VB*4 ,27*( R21**2 . -R1 1**2 . ) 

VW a 4 .27* ( R22**2 . -R12**2 . ) 

WTCMC»RH01*VB+RH02*VW 
WRITE(11,*)' 

WRITE(11,*) 

WRITE(11,*)' 

WRITEdl,*)' 

WRITE(11,*)' 

WRITEdl,*) ' 

WRITE(11,*) 

WRITE(11,*) 

WRITE(11,*) 

GO TO 10 
READ(10,*)D 
READ(10,*)RH0P 
READ(10,*)S 
READ(10,*)RH0 
WRITE(11,*)' 

WRITE(11,*) 

WRITE(11,*)' 

WRITE(11,*) 

WRITE(11,*) ' 


OUTPUT' 

BUMPER THICKNESS = ' ,T1 , ' CM ' 
WALL THICKNESS = ' ,T2 , ' CM* 

MIN. WEIGHT a ' ,WT , ' CM* 

CMC MIN. WEIGHT a ',WTCMC,'KG' 


MADDEN' 

INPUT' 

PROJECTILE DIAMETER IN CM = ' ,D 
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0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 46 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 55 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 


WRITE (11 ,*) 1 PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 

WRITE( 11 ,*) ' BUMPER/WALL DENSITY IN GM/CUBIC CM = ' ,RHO 

WRITE (11 ,*) 1 BUMPER/WALL SEPARATION IN CM = ' ,S 

WRITE (11 ,*) 

T1T=0.0 
T2T=0.0 
DO 46 J=1 ,16 
V=FLOAT( J ) 

CALL MADDEN ( V , D , RHOP , S ,RHO ,T1 , T2 , WT , WTCMC ) 

T1T=T1T+T1*XPV(J) 

T2T=T2T+T2*XPV(J ) 

CONTINUE 
T1=T1T 
T2=T2T 
WT=T1+T2 
R12=211 . 

R22=211.+T2 

Rll=211 .+T2+S 

R21=211 .+T1+T2+S 

VB=4 ,27*( R2 1**2. -R1 1**2.) 

VW=4 .27* (R22**2.-R 12**2 . ) 

WTCMC=RHO*( VB+VW) 

WRITE(11,*)' OUTPUT' 

WRITE(11,*) 

WRITE (11 ,*) ' BUMPER THICKNESS = ' ,T1 , ' CM ' 

WRITE (11 ,*) 1 WALL THICKNESS = , ,T2, , CM' 

WRITE(11 ,*) 1 MIN. WEIGHT = ',WT,'CM' 

WRITE (11 ,*) 1 CMC MIN. WEIGHT = '.WTCMC, 'KG' 

WRITEdl ,*) 

WRITE (11,*) 

WRITE(11,*) 

GO TO 10 

READ( 10 ,*)D 

READ(10,*)RH0P 

READ(10,*)RH01 

READ(10,*)RH02 

READ(10 ,*) S 

READ(10,*)XL2 

WRITE(11,*)' WILKINSON' 

WRITE(11 ,*) 

WRITE(11,*)' INPUT' 

WRITE(11,*) 

WRITE(11,*)' PROJECTILE DIAMETER IN CM = ' ,D 

WRITE(11,*)' PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 

WRITEdl,*)' BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 

WRITE(11,*)' WALL DENSITY IN GM/CUBIC CM = ' ,RH02 

WRITE (11,*) ' BUMPER/WALL SEPARATION IN CM = ' ,S 

WRITE(11,*)' WALL MATERIAL CONSTANT = ' ,XL2 

WRITE(11,*) 

T1T=0.0 
T2T=0.0 
DO 56 J=1 ,16 
V=FLOAT( J ) 

CALL W I LK I NSON ( V ,D ,RH0P,RH01,RH02,S,XL2 , 

T1,T2,WT, WTCMC) 

T1T=T1T+T1*XPV( J) 

T2T=T2T+T2*XPV(J ) 
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0172 

56 


CONTINUE 

1 

0173 



T1=T1T 

m 

0174 



T2=T2T 


0175 



WT S T1+T2 

1 

0176 



R12=211 . 


0177 



R22=211.+T2 

i 

0178 



R11-211.+T2+S 


0179 



R21=211 .+T1+T2+S 

i 

0180 



VB=4.27*(R21**2.-R11**2.) 

0181 



VW=4.27*(R22**2.-R12**2 .) 

0182 



WTCMC=RH01*VB+RH02 

*vw 1 

0183 



WRITE(11,*)' 

OUTPUT' 

0184 



WRITE (11,*) 

1 

0185 



WRITE (11 ,*) ' 

BUMPER THICKNESS = ' , T1 , ' CM ' 

0186 



WRITE (11 ,*) 1 

WALL THICKNESS - ' ,T2,'CM' 1 

0187 



WRITE(11,*)' 

MIN. WEIGHT = ' ,WT,'CM' 

0188 



WRITE (11 ,*) 1 

CMC MIN. WEIGHT = ',WTCMC,'KG' j 

0189 



WRITE(11,*) 

1 

0190 



WRITE (11,*) 


0191 



WRITE (11,*) 

I 

0192 



GO TO 10 


0193 

65 


READ(10,*)D 

1 

0194 



READ(10 ,*)RH01 


0195 



READ(10,*)S 


0196 



READ(10,*)THETA 

• 

0197 



READ(10,*)XN 


0198 



READ (10,*) El 

i 

0199 

★ ★*** 

MODIFIED BURCH ***** 

0200 



WR ITE( 11 ,*) ' 

MODIFIED BURCH' 1 

0201 



WRITE( 11 ,*) 


0202 



WRITE(11 ,*) ' 

INPUT' 

0203 



WRITE(11 ,*) 


0204 



WRITE( 11 ,*) 1 

'PROJECTILE DIAMETER IN CM = ' ,D 

0205 



WRITE( 11 ,*) 1 

BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 

0206 



WR I TE (11,*) 1 

BUMPER/WALL SEPARATION IN CM * ' ,S 

0207 



WRITE (11 ,*) ' 

IMPACT ANGLE FROM NORMAL IN DEGREES = '.THETA 

0208 



WRITE(11 ,*) 1 

NUMBER OF PLATES TO PENETRATE AFTER FIRST' , 

0209 


& 

1 

BUMPER = ' ,XN 

0210 



WRITE(11,*)' 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE ' , 

0211 


& 


' SEC = ' ,E1 

0212 



WRITE(11,*) 


0213 



T1T=0.0 


0214 



T2T=0.0 


0215 



DO 66 J=1 ,16 


0216 



V=FLOAT( J ) 

1 

0217 



CALL BURCH(V,D,RH01,S, THETA, 1 

0218 


& 

XN,E1 , 

T1,T2,WT,WTCMC) 

0219 



T1T=T1T+T1*XPV( J) 


0220 



T2T=T2T+T2*XPV( J ) 


0221 

66 


CONTINUE 


0222 



T1=T1T 


0223 



T2=T2T 


0224 



WT=T1+T2 


0225 



R12=211 . 


0226 



R22=211.+T2 


0227 



Rll=211 .+T2+S 


0228 



R21=211 .+T1+T2+S 
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0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 10 

0243 

0244 


VB=4. 27* (R2 1**2 .-Rl 1**2 .) 
VW=4.27*(R22**2.-R12**2.) 

WTCMC=RH01*VB+2.81*VW 
WRITE(11 ,*) ' OUTPUT' 

WRITE (11 ,*) 

WR I TE (11 [*) ' BUMPER THICKNESS = ' ,T1 , ' CM ' 

WRITE (11,*) ' WALL THICKNESS = , ,T2, , CM I 

WRITE( 11 ,*) ' MIN. WEIGHT = ' ,WT, ' CM' 

WRITE(11 ,*) 1 CMC MIN. WEIGHT = '.WTCMC.'KG' 

WRITE(11 ,*) 

WRITE (11 ,*) 

WRITE(11 ,*) 

GO TO 10 
CONTINUE 
STOP 
END 


PROGRAM SECTIONS 



Name 

Bytes 

Attributes 







0 

$CODE 

4728 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

1 

SPDATA 

741 

PIC CON REL 

LCL 

SHR 

NOEXE 

RD 

NOWRT 

LONG 

2 

SLOCAL 

692 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RO 

WRT 

LONG 


Total Space Allocated 

6161 









ENTRY POINTS 

Address Type Name 
0-00000000 IMPACT5V$MAIN 


VARIABLES 


Address 

Type 

Name 

Address 

Type 

Name 

Address 

Type 

Name 

2-0000004C 

R*4 

D 

2-0000008C 

R*4 

El 

2-00000050 

R*4 

H 

2-00000040 

1*4 

IV 

** 

1*4 

J 

2-00000048 

1*4 

NCODE 

** 

R*4 

Rll 

** 

R*4 

R12 

** 

R*4 

R21 

2-00000090 

R*4 

RHO 

2-0000006C 

R*4 

RH01 

2-00000070 

R*4 

RH02 

2-00000074 

R*4 

S 

2-0000007C 

R*4 

SY1 

2-00000080 

R*4 

SY2 

** 

R*4 

TIT 

2-0000005C 

R*4 

T2 

** 

R*4 

T2T 

2-00000054 

R*4 

V 

** 

R*4 

VB 

** 

R*4 

VW 

2-00000064 

R*4 

WTCMC 

2-00000078 

R*4 

XL2 

2-00000088 

R*4 

XN 
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0001 SUBROUTINE NYSMITH( V,D,H,T1 ,T2,WT,WTCMC) 

0002 DMAX=0.24*H*V**-0.2 

0003 IF(D.GT.DMAX)THEN 

0004 WRITE(11 ,*) ' NO SOLUTION— PROJ. DIA. TOO LARGE FOR NYSMITH' 

0005 GO TO 16 

0006 ENDIF 

0007 T1=1.93*V**0.18*D**1.91/H**0.91 

0008 T2=1.86*T1 

0009 16 CONTINUE 

0010 RETURN 

0011 END 


VAX 

SAI 


FORTR 

USRD 






PROGRAM SECTIONS 



Name 

Bytes 

Attributes 







0 

$CODE 

150 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

1 

$PDATA 

49 

PIC CON REL 

LCL 

SHR 

NOEXE 

RD 

NOWRT 

LONG 

2 

$LOCAL 

8 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 


Total Space Allocated 

207 









ENTRY POINTS 


Address Type 

Name 






0-00000000 

NYSMITH 






VARIABLES 







Address Type 

Name 

Address 

Type 

Name 

Address Type 

Name 

AP-000000080 R*4 

D 

** 

R*4 

DMAX 

AP-OOOOOOOCO R*4 

H 

AP-000000140 R*4 

T2 

AP-00000004@ 

R*4 

V 

AP-000000180 R*4 

WT 


LABELS 

Address Label 
0-00000095 16 
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cat ncnnrc 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


BOEING 


iu-nuy- 130 / n.tt.ju om 

SUBROUTINE B0EI NG( V ,D ,RH0P ,RH01 ,RH02 ,S ,XL2 ,SY1 ,SY2 .THETA , 

& XN.El ,T1 .T2.WT.WTCMC) 

***** PEN4 ***** 

Tl=0.16 
V=V* 3280. 

D=D/30.48 
RP=D/2 .0 
RH0P=RH0P*1 .94 
RH01=RH01*1 .94 
RH02=RH02*1 .94 
NITSP=0 
NITSP=NITSP+1 
NP1=0 

TlP s Tl/30.48 

T2P=FT2P(RH0P ,V ,RP ,SY1 .THETA.RH01 ,SY2 ,D ,RH02 ,T1P) 

WT=T1P+T2P 

IF(NITSP.EQ.1)THEN 

T1P1=1.1*T1P 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 

WT1*T1P1+T2P1 

ENDIF 

IF(WT1.GT.WT)THEN 

T1P1*0.82*T1P1 

T2P1=FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 

WT1=T1P1+T2P1 

590 IF(WT1.GT.WT)THEN 

GO TO 601 
ELSE 

T1P-T1P1 

T2P a T2Pl 

WT=WT1 

NP1=NP1+1 

IF(NP1.EQ.100)THEN 

WRITE(11 ,*) ' NO CONVERGENCE IN PEN4 1 
GO TO 557 
ENDIF 

TlPl a 0.9*TlPl 

T2P1=FT2P(RH0P,V,RP,SY1 .THETA, RH01 ,SY2,D ,RH02, TIPI ) 

WT1*T1P1+T2P1 
GO TO 590 
ENDIF 
ELSE 

579 T1P=T1P1 

T2P=T2P1 
WT=WT1 

T1P1 S 1 .1*T1P1 

T2P1=FT2P(RH0P,V ,RP,SY1, THETA, RH01 ,SY2,D,RH02, TIPI) 

WT1=T1P1+T2P1 
IF(WT1.GT.WT)THEN 
GO TO 601 
ELSE 

NP1=NP1+1 

I F ( NP1 .EQ.100)THEN 

WRITE (11 ,*) ' NO CONVERGENCE IN PEN4 1 
GO TO 557 
ENDIF 
GO TO 579 


ORIGINAL PAGE IS 
QE POOR QUALITY 
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uorwj i of 


FORTRAN 

i icno t ci/ 


iO-Muy-iso/ u.tc.au 


OMi uor\u 


0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 

0113 

0114 


601 


557 

***** 


***** 


125 


1099 


ENDIF 
ENDIF 
CONTINUE 
0=30. 48*D 
RH0P=RH0P/1 .94 
RH01=RH0l/l .94 
RH02=RH02/1 .94 
T1P=30.48*T1P 
T2P=30.48*T2P 
IF(T1P/D.LE.0.4)VF=4100 
IF(T1P/D.GT.0.4)VF=4986*(T1P/D)**0.21 


VF=VF+4000. 

IF(V.LE.VF)THEN 

WRITE(11,*)' INSIDE OF PEN4 LIMITS' 
T1=T1P 
T2=T2P 
GO TO 1102 
ENDIF 
CONTINUE 

WILKINSON ***** 

V=V/3280. 

T1=0.604*D**2.*RH0P/S 
T1=T1*SQRT(V/(XL2*RH01*RH02) ) 
RATI0=D*RH0P/(T1*RH01) 

IF( RATIO. GT.1.0)T2=T1 
IF(RATI0.LE.1.0)T2=T1/RATI0 
MODIFIED BURCH ***** 

V B=V*3280. 

DB=D/2.54 

CM=SQRT(E1/RH01) 

CM=CM/30.48 

SB=S/2.54 

IF(TH£TA,LE.0.001)G0 TO 125 
CHI=TAN( THETA) -0.5 

F2=0.5-1.87*(T1B/D)+(5.*T1B/D-1.6)*CHI**3.0 

F2=F2+(1.7-12.*T1B/D)*CHI 

F3=0.32*(T1B/D)**0.83 

F3=F3+0.48*(T1B/D)**0.33*(SIN(THETA) )**3.0 

T2F=D* ( ( Fl+O .63*F2 ) /XN ) *( CM/ V ) **2 .29 

T2F=T2F*(D/S)**0 .71 

T2N=F3*( CM/ V ) **1 .33*D/XN 

IF(T2N.GE.T2F)T2B=T2N 

IF(T2N.LT.T2F)T2B=T2F 

T2B=T2B*2 .54 

IF(T2B.GT.T2)NREGI0N=3 

IF(T2B.GT.T2)T2=T2B 

GO TO 155 

CONTINUE 

TlBl=Tl/2.54 

NITSB=0 

XK1=(DB/XN)**1.71*(CM/VB)**2.29/SB**0.71 


original page rs 

OE poor quality 


VDELTA=0.0 

DELTA3=0.52 

DELTA2=2.33*(1.-1.57*DELTA3) 

DELTA1=1 .33*(2 .*DELTA3-1 .) 

VDELTA1=( 1 ./DELTA1 ) **DELTA1*( 2 ,8*XK1/ ( DELTA2*DB**0 .57 ) )**DELTA2 
VDELTA1=VDELTA1*(1.58*XK1*DB**0.57/DELTA3)**DELTA3 


10-Aug-1987 11:43:04 

in a . . ^ inoi ii.Ao.cr 


VAX FORT 




cat itcnni c 


BOEING 


iu-Huy-j.»o/ u;tt;3u 


oh i uoruoior 


0115 


IF(VDELTA1.LT.VDELTA)THEN 

0116 


DELTA1=1 .33*(2.*DELTA3-1.04) 

0117 


T1B=DELTA1*VDELTA 

0118 


T2B*VDELTA-T1B 

0119 


GO TO 499 

0120 


ENDIF 

0121 


VDELTA=VDELTA1 

0122 


DELTA3=DELTA3+0.02 

0123 


IF(DELTA3.GT.0.63)THEN 

0124 


T1B=DELTA1*VDELTA 

0125 


T2B=VDELTA-T1B 

0126 


GO TO 499 

0127 


ENDIF 

0128 


GO TO 1099 

0129 

499 

CONTINUE 

0130 

***** 

COMPARISON OF MODIFIED BURCH AND WILKINSON 

0131 

199 

CONTINUE 

0132 


T10W=Tl/2.54 

0133 


F10W=1.58*(DB/T10W)**0.57+2.80*(T10W/DB)**0 

0134 


T2BT10W=(F10W/XN)**1.71*(CM/VB)**2.29*DB**1 

0135 


T2BT10W=T2BT10W/SB**0 .71 

0136 


T2BT10W=T2BT10W*2.54 

0137 


RATI0B=(DB*RH0P)/(RH01*T1B) 

0138 


T2WT10B=0 ,364*D**3 .*RH0P*V/ ( XL2*RH02*S**2 . ) 

0139 


IF(RATI0B .GT .1 .0 )T2WT10B=T2WT10B*RATI0B 

0140 


IF(T2BT10W.GT.T2)T2=T2BT10W 

0141 


T2B=T2B*2.54 

0142 


I F( T2WT1 OB .GT .T2B ) T2B=T2WT10B 

0143 


T1B=T1B*2.54 

0144 


IF(T1B+T2B.LT.T1+T2)THEN 

0145 


T1=T1B 

0146 


T2=T2B 

0147 


ENDIF 

0148 

155 

CONTINUE 

0149 

1102 

IF(T2.LE.0.01)THEN 

0150 


SIGMA=SY2/144. 

0151 


T2=3099. 1/SIGMA 

0152 


ENDIF 

0153 

156 

RETURN 

0154 


END 


***** 
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in An** mm 11.A0.ct cat 


FORTRAN 

ucnn r ct/ 


BOEING 



xu-Muy-i^u/ n.^c.uu 


onx uorwj; 


0001 FUNCTION FT2B(DB,T1B,XN,CM,VB,SB) 

0002 F1=2.42*(DB/T1B)**0.33+4.26*(T1B/DB)**0 .33 

0003 Fl-Fl-4.18 

0004 FT2B-(F1/XN)**1.71*(CM/VB)**2.29*DB**1.71/SB**0.71 

0005 RETURN 

0006 END 


PROGRAM SECTIONS 


Name 

0 $C0DE 
2 $L0CAL 

Total Space Allocated 


Bytes 

Attributes 







163 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

4 

167 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 


ENTRY POINTS 


Address Type Name 
0-00000000 R*4 FT2B 

VARIABLES 


Address Type 

Name 

Address Type 

Name 

Address 

Type 

Name 

AP -000000 100 R*4 

CM 

AP -000000040 R*4 

DB 

** 

R*4 

FI 

AP -000000080 R*4 

TIB 

AP-000000140 R*4 

VB 

AP-OOOOOOOC0 

R*4 

XN 


10-Aug-1987 11:43:04 

lO A..^ 1AOT 11 .AO.CC 


VAX 

CAT 


FORTH 

1 1 cnn t c 


iu-ttuy-i 30/ 


om uorwi-ir 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 999 

0016 


FUNCTION FT2P(RH0P ,V ,RP ,SY1 , THETA, RH01 ,SY2 ,D,RH02 .TIP) 

A=1.33*RH0P*(V*RP)**2. 

B=8.0*SY1*EXP(-3.125E-04*V)/C0S(THETA) 

C=1.33*RH0P*RP**2.0 
D1=RP*RH01/C0S( THETA) 

XK1«1.67*(RH0P/(2.*SY2))**0.31 
XK1=XK1*(0 .281 *D*RH0P/RH02) **0.33 
XK1=XK1*C0S( THETA) 

C1P1 S (A-B*T1P)/(C+D1*T1P) 

IF(C1P1.LE. 0.001) THEN ORIGIN ' L 

FT2P-0.0 
GO TO 999 
ENDIF 

FT2P=XK1*C1P1**0 .31 
RETURN 
END 




OE POOR QUALITY 


PROGRAM SECTIONS 


Name 

Bytes 

Attributes 







0 $CODE 

222 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

2 $LOCAL 

4 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 

Total Space Allocated 

226 









ENTRY POINTS 

Address Type Name 
0-00000000 R*4 FT2P 


VARIABLES 


Address 

Type 

Name 

Address 

Type 

Name 

Address Type 

Name 

** 

R*4 

A 

★* 

R*4 

B 

** R*4 

C 

AP-00000020@ 

R*4 

D 

** 

R*4 

D1 

AP-00000018@ R*4 

RH01 

AP-00000004@ 

R*4 

RHOP 

AP-0000000C@ 

R*4 

RP 

AP-00000010® R*4 

SY1 

AP-00000028@ 

R*4 

TIP 

AP-00000014@ 

R*4 

THETA 

AP-00000008@ R*4 

V 


LABELS 

Address Label 
0-000000D9 999 


FT2P 


10-Aug-1987 11:43:04 VAX FORTRAN 

in i no-/ ii.4o.cc cat iicontci/ 


lU-Muy- 130 / 


OMl U3t\V\ 


0001 SUBROUTINE MADDEN(V,D,RH0P,S,RH0,T1 ,T2,WT,WTCMC) 

0002 V=V*1.E05 

0003 T1=0 ,009*SQRT( V ) *RH0P*D**2 .0 

0004 Tl=Tl/( S*RH0**1 .5) 

0005 T2=T1 

0006 RETURN 

0007 END 


PROGRAM SECTIONS 


Name 

Bytes 

Attributes 

0 SCODE 

69 

PIC CON REL LCL SHR EXE RD NOWRT LONG 

Total Space Allocated 

69 



ENTRY POINTS 


Address Type 

Name 





0-00000000 

MADDEN 





VARIABLES 






Address Type 

Name 

Address Type 

Name 

Address Type 

Name 

AP-00000008@ R*4 

D 

AP-00000014@ R*4 

RHO 

AP-OOOOOOOC0 R*4 

RHOP 

AP-00000018@ R*4 

T1 

AP-0000001C@ R*4 

T2 

AP-00000004@ R*4 

V 

AP-00000024@ R*4 

WTCMC 






FUNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTHSSQRT 


10-Aug-1987 11:43:04 VAX FORTW 

in a..~ i noi ii.ao.cc cat ncr>nt< 


J.u-«uy- 130 / li.Hd.JU 


on i uorujior 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 
0009 


PROGRAM SECTIONS 
Name 
0 $C0DE 

Total Space Allocated 


ENTRY POINTS 

Address Type Name 
0-00000000 WILKINSON 


VARIABLES 


Address Type 

Name 

Address Type 

Name 

Address Type 

Name 

AP -000000080 R*4 

D 

** R*4 

RATIO 

AP-000000100 R*4 

RH01 

AP-OOOOOOOC0 R*4 

RHOP 

AP-000000180 R*4 

S 

AP-000000200 R*4 

T1 

AP -000000040 R*4 

V 

AP-000000280 *R*4 

WT 

AP-OOOOOO2C0 R*4 

WTCMC 


FUNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTHSSQRT 


Bytes Attributes 

87 PIC CON REL LCL SHR EXE RD NOWRT LONG 
87 


SUBROUTINE WILKINS0N(V,D,RH0P,RH01,RH02,S,XL2, 
T1 ,T2 ,WT,WTCMC) 
T1=0.604*D**2.*RH0P/S 
T1=T1*SQRT(V/(XL2*RH01*RH02) ) 
RATI0=D*RH0P/(T1*RH01) 

IF(RATI0.GT.1.0)T2=T1 

IF(RATI0.LE.1.0)T2=T1/RATI0 

RETURN 

END 1 
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0001 


MODIFIED BURCH ***** 

0002 


SUBROUTINE BURCH( V ,D ,RH01 ,S, THETA, 

0003 

& 

XN,E1 ,T1 ,T2,WT,WTCMC) 

0004 


VB=V*3280. 

0005 


DB=D/2 .54 

0006 


CM=SQRT(E1/RH01) 

•0007 


CM=CM/30.48 

0008 


SB=S/2 .54 

0009 


IF(THETA.LE.0.001)G0 TO 425 

0010 


CHI=TAN( THETA) -0.5 

0011 


F2=0.5-1.87*(T1B/D)+(5.*T1B/D-1.6)*CHI**3.0 

0012 


F2=F2+(1.7-12.*T1B/0)*CHI 

0013 


F3=0.32*(T1B/D)**0 .83 

0014 


F3=F3+0.48*(T1B/D)**0 .33*( SIN(THETA) )**3 .0 

0015 


T2F=D*( ( Fl+0 ,63*F2)/XN)*(CM/V)**2 .29 

0016 


T2F*T2F*(D/S)**0.71 

0017 


T2N=*F3*(CM/V)**1 ,33*D/XN 

0018 


IF{T2N.GE.T2F)T2B=T2N 

0019 


IF(T2N.LT.T2F)T2B=T2F 

0020 


T2B=T2B*2 .54 

0021 


IF(T2B.GT.T2)NREGI0N=3 

0022 


IF(T2B.GT.T2)T2=T2B 

0023 


GO TO 499 

0024 

425 

CONTINUE 

0025 


NITSBO 

0026 


XK1=(DB/XN)**1 .71* (CM/VB) **2 .29/ SB**0 .71 

0027 


VDELTA=0.0 

0028 


D£LTA3=0 .52 

0029 

1099 

DELTA2=2 .33*( 1 .-1 ,57*DELTA3 ) 

0030 


DELTA1=1 .33*(2.*DELTA3-1. ) 

0031 


VDELTA1=( 1 ./DELTA1 )**DELTA1*(2.8*XK1/(DELTA2*DB**0.57 ) )**DELTA2 

0032 


VDELTA1=VDELTA1*(1.58*XK1*DB**0.57/DELTA3)**DELTA3 

0033 


IF(VDELTA1.LT.VDELTA)THEN 

0034 


DELTA1=1.33*(2 ,*DELTA3-1 .04) 

0035 


T1=DELTA1*VDELTA 

0036 


T2=VDELTA-T1 

0037 


GO TO 499 

0038 


ENDIF 

0039 


VDELTA=VDELTA1 

0040 


DELTA3=DELTA3+0 .02 

0041 


IF(DELTA3.GT.0.63)THEN 

0042 


T1=0ELTA1*VDELTA 

0043 


T2 S VDELTA-T1 

0044 


GO TO 499 

0045 


ENDIF 

0046 


GO TO 1099 CK1GMAL PAGE IS 

0047 

0048 

499 

Tl-Tl‘2^54 

0049 


T2=T2*2.54 

0050 


RETURN 

0051 


END 
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IMPACT5VM 


IMPACT5VM is a spacecraft protective systems design optimization code 
similar to IMPACT5V. IMPACT5VM differs from IMPACT5V in that the 
optimal design is weighted according to the chosen meteoroid velocity 
probability distribution. IMPACT5VM employs the Geometric 
Programming optimization technique in evaluating a number of nonlinear 
piecewise continuous, functional impact predictors. These include the 

Nysmith, Boeing, Madden, Wilkinson, and Burch predictors. Inputs vary 
depending on the predictor used, however, typical inputs include the 
meteoroid velocity distribution file, projectile characteristics (as 

determined from METEOR 1 or DEBRIS 1), design material properties, and 
general design configuration. In particular, IMPACT5VM is an 
optimization code for a single bumper, single wall configuration. Outputs 
include the optimal design thicknesses (bumper and wall) and the 

minimum design weight. The velocity distribution file for meteoroids, 
(METVEL.IN), follows. The file being used must be assigned to FOR012 
before running IMPACT5VM. Sample input (IMPACT5VM.IN), output 
(IMPACT5VM.OUT), and program listing (IMPACT5VM.LIS) follow the 

velocity distribution file. 





9 

.045 

10 

.045 

11 

.045 

12 

.045 

13 

.25 

14 

.25 

15 

.25 

16 

.26 

17 

.26 

18 

.26 

19 

.14 

20 

.14 

21 

.14 

22 

.10 

23 

.10 

24 

.10 

2b 

.06 

26 

.06 

27 

.06 

28 

.05 

29 

.05 

30 

.05 

31 

.025 

32 

.025 

33 

.025 

34 

.015 

35 

.015 

36 

.015 

37 

.012 

38 

.012 

39 

.012 

40 

.011 

41 

.011 

42 

.011 

43 

.009 

44 

.009 

45 

.009 

46 

.0075 

47 

.0075 

48 

.0075 

49 

.0075 

50 

.0075 

51 

.0075 

52 

.0075 

53 

.0075 

54 

.0075 

55 

.0075 

56 

.0075 

57 

.0075 

58 

.0075 

59 

.0075 

60 

.0075 

61 

.0075 

62 

.0075 

63 

.0075 

64 

.0075 

65 

.0075 

66 

.0075 

67 

.0075 

68 

.0075 



69 

70 

71 

72 


.0075 

.0075 

.0075 

.0075 





5 

1 

0.84 

10 . 

2 

0.84 

2.81 

2.81 

2.81 

10 . 

0.401 

7344000. 

7344000. 

0 . 

0.85 

7.239E11 

3 

0.84 

2.81 

10 . 

2.81 

4 

0.84 

2.81 

2.81 

2.81 

10 . 

0.401 

5 

0.84 

2.81 

10 . 

0 . 

0.85 

7.239E11 


NUMBER OF CASES 

NYSMITH PREDICTOR 

PROJ. DIAMETER IN CM 

BUMPER/WALL SEPARATION IN CM 

BOEING PREDICTOR 

PROJ. DIA. IN CM 

PROJ. DENSITY IN GM/CUBIC CM 

BUMPER DENSITY IN GM/CUBIC CM 

WALL DENSITY IN GM/CUBIC CM 

BUMPER/WALL SEPARATION IN CM 

WALL MATERIAL CONSTANT 

BUMPER YIELD STRENGTH IN LB/SQUARE FOOT 

WALL YIELD STRENGTH IN LB/SQUARE FOOT 

IMPACT ANGLE FROM NORMAL 

NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 

YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 

MADDEN PREDICTOR 

PROJ. DIA. IN CM 

PROJ. DENSITY IN GM/CUBIC CM 

BUMPER/WALL SEPARATION IN CM 

BUMPER/WALL DENSITY IN GM/CUBIC CM 

WILKINSON PREDICTOR 

PROJ. DIA. IN CM 

PROJ. DENSITY IN GM/CUBIC CM 

BUMPER DENSITY IN GM/CUBIC CM 

WALL DENSITY IN GM/CUBIC CM 

BUMPER/WALL SEPARATION IN CM 

WALL MATERIAL CONSTANT 

BURCH PREDICTOR 

PROJ. DIA. IN CM 

BUMPER DENSITY IN GM/CUBIC CM 

BUMPER/WALL SEPARATION IN CM 

IMPACT ANGLE FROM NORMAL 

NUMBER OF PLATES TO PENETRATE AFTER 1ST BUMPER 
YOUNG'S MODULUS FOR BUMPER IN GM/CM-SQUARE SEC 




NYSMITH 


INPUT 

PROJECTILE DIAMETER IN CM * 0.8400000 

BUMPER/WALL SEPARATION IN CM * 10.00000 

OUTPUT 

BUMPER THICKNESS * 0.2964782 CM 

WALL THICKNESS * 0.5514496 CM 

MIN. WEIGHT = 0.8479279 CM 

CMC MIN. WEIGHT * 4373.689 KG 


BOEING 

INPUT 

PROJECTILE DIAMETER IN CM - 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

BUMPER YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

WALL YIELD STRENGTH IN LB/SQUARE FT = 7344000. 

IMPACT ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 
NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7 .2389997E+11 

OUTPUT 

BUMPER THICKNESS = 0.3065161 

WALL THICKNESS = 0.3286171 

MIN. WEIGHT = 0.6351332 CM 

CMC MIN. WEIGHT = 3294.346 


MADDEN 

INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

OUTPUT 

BUMPER THICKNESS = 0.5543199 CM 

WALL THICKNESS * 0.5543199 CM 

MIN. WEIGHT = 1.108640 CM 

CMC MIN. WEIGHT = 5761.344 KG 


CM 

CM 

KG 


WILKINSON 



INPUT 


PROJECTILE DIAMETER IN CM = 0.8400000 

PROJECTILE DENSITY IN 6M/CUBIC CM = 2.810000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

WALL DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

WALL MATERIAL CONSTANT = 0.4010000 

OUTPUT 

BUMPER THICKNESS = 0.3114124 CM 

WALL THICKNESS = 0.3114124 CM 

MIN. WEIGHT » 0.6228247 CM 

CMC MIN. WEIGHT = 3233.087 KG 


MODIFIED BURCH 
INPUT 

PROJECTILE DIAMETER IN CM = 0.8400000 

BUMPER DENSITY IN GM/CUBIC CM = 2.810000 

BUMPER/WALL SEPARATION IN CM = 10.00000 

IMPACT ANGLE FROM NORMAL IN DEGREES = O.OOOOOOOE+OO 

NUMBER OF PLATES TO PENETRATE AFTER FIRST BUMPER = 0.8500000 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE SEC = 7 .2389997E+11 

OUTPUT 

BUMPER THICKNESS * 4.5066189E-02CM 

WALL THICKNESS = 9.7935580E-02CM 

MIN. WEIGHT = 0.1430018 CM 

CMC MIN. WEIGHT = 735.0610 KG 


1 
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17-Aug-1987 15:19:56 

SAI_USRDISK 

1 )001 


DIMENSION XPV(IOO) 


1 >002 


UPEN(UNIT*12,TYPE='0LD' .ACCESS® 'SEQUENTIAL' ) 

OPEN( UNIT=10 .TYPE* ' OLD ' .NAME* ' IMPACT5VM.IN ' .ACCESS* 1 SEQUENTIAL ' ) 


P )003 



J004 


OPEN ( UNI T*ll .TYPE* 'NEW .NAME* ' IMPACT5VM.0UT' .ACCESS*' SEQUENTIAL ' ) 


| J005 


DO 24 1*1,64 


1 )006 


READ(12,*) IV.XPV(IV) 


J007 

24 

CONTINUE 


1 )008 


READ( 10 ,* ) NRUNS 


I )009 


DO 10 1=1, NRUNS 


1 )01U 


READ(10,*)NC0DE 


)011 


IF(NC0DE.EQ.1)G0 TO 25 


| )012 


I F( NCODE .EQ.2)G0 TO 35 


1 JO 1 3 


IF(NC0DE.EQ.3)G0 TO 45 


)014 


I F( NCOOE .EQ.4)G0 TO 55 


■ JO 1 5 


IF(NC0DE.EQ.5)G0 TO 65 


I JO 1 6 

25 

READ(10,*)D 


1 )017 


READ(10,*)H 


J018 


WRITE(11,*)' NYSMITH' 


J019 


WRITEdl,*) 


J020 


WRITEdl,*)' INPUT' 


J021 


WRITEdl,*) 


. J022 


WRITEdl,*)' PROJECTILE DIAMETER IN CM = ' ,D 


)023 


WRITE(11,*)' BUMPER/WALL SEPARATION IN CM * ' ,H 


1 J024 


WRITE(11,*) 


)02b 


T1T=0.0 


1 J026 


T2T=U.O 


|- )027 


DO 26 J=1 ,64 


)028 


K=J+8 


)029 


V*FLOAT(K) 


)030 


CALL NYSMITH(V,D,H,T1 ,T2,WT,WTCMC) 


)031 


T1T=T1T+T1*XPV(K) 


)032 


T2T=T2T+T2*XPV(K) 


1 )033 

26 

CONTINUE 


| )034 


T1=T1T/3.1335 


)035 


T2=T2T/3.1335 


, )036 


WT*T1+T2 


>037 


WTCMC=Tl**2.+2 .*T1*(211.+T2+H)+T2**2.+422.*T2 


l )038 


WTCMC*12.*WTCMC 


)039 


WRITE(11,*)' OUTPUT' 


| )040 


WRITE (11 ,*) 


>041 


WRITE(11,*)' BUMPER THICKNESS * '.Tl.'CM' 


)042 


WRITEdl,*)' WALL THICKNESS = ',T2,'CM' 


)043 


WRITE(11,*)' MIN. WEIGHT = ',WT,'CM' 


)U44 


WRITEdl,*)' CMC MIN. WEIGHT = '.WTCMC.'KG' 


)045 


WRITE(11,*) 


)046 


WRITE(11,*) 


>047 


WRITE(11,*) 


)048 


GO TO 10 


>049 

35 

READ(10,*)D 


>050 


READ(10,*)RH0P 


)051 


READ(10,*)RH01 


)052 


READ(10,*)RH02 


)0b3 


READ(10,*)S 


)054 


READ(10,*)XL2 


)055 


READ(10,*)SY1 


J056 


READ(10,*)SY2 


| J057 

1 


READ(10,*)THETA 



[MPACT5VMSMAIN 
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3058 



READ(10,*)XN 


3059 



READ(10,*)E1 


)060 



WRITE (11 ,*) 1 

BOEING' 

J06 1 



WRITE(11 ,*) 


)U62 



WRITE (11 ,*) ' 

INPUT' 

)063 



WRITE(11 ,*) 


)064 



WRITEdl,*)' 

PROJECTILE DIAMETER IN CM 3 ' ,D 

)065 



WRITEdl,*)' 

PROJECTILE DENSITY IN GM/CUBIC CM 3 ' ,RHOP 

)066 



WRITE(11 ,*) ' 

BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 

)U67 



WRITEdl,*)' 

WALL DENSITY IN GM/CUBIC CM = ' ,RH02 

)068 



WRITEdl,*)' 

BUMPER/WALL SEPARATION IN CM = ' ,S 

3069 



WRITE(11,*)' 

WALL MATERIAL CONSTANT 3 ' ,XL2 

)070 



WRITEdl ,*) ' 

BUMPER YIELD STRENGTH IN LB/SQUARE FT 3 ' ,SY1 

3071 



WRITE(11,*)' 

WALL YIELD STRENGTH IN LB/SQUARE FT = ' ,SY2 

3072 



WRITE(11,*) ' 

IMPACT ANGLE FROM NORMAL IN DEGREES 3 '.THETA 

)073 



WRITEdl,*)' 

NUMBER OF PLATES TO PENETRATE AFTER FIRST', 

)074 


& 


' BUMPER 3 ' ,XN 

3075 



WRITEdl,*)' 

BUMPER YOUNGS MODULUS IN GM/CM-SQUARE' , 

3076 


& 


' SEC 3 ' ,E1 

)077 



WRITE (11 ,*) 


3078 



T1T 3 0.0 


3079 



T2T 3 0.0 


)080 



DO 36 J=1 ,64 


0081 



K=J+8 


3082 



V=FLOAT(K) 


)083 



CALL B0EING(V,D,RH0P,RH01,RH02,S,XL2,SY1,SY2, THETA, 

)084 


& 

XN,E1 

,T1,T2,WT,WTCMC) 

)085 



T1T=T1T+XPV ( K) *T1 


3086 



T2T=T2T+XPV(K)*T2 


)087 

36 


CONTINUE 


)088 



T1=T1T/3.1335 


)089 



T2=T2T/3.1335 


)090 



WT=T1+T2 


)091 



R12*211 . 


JU92 



R22=211.+T2 


)093 



Rll 3 211 .+T2+S 


)(J94 



R21®211 .+T1+T2+S 


)09b 



VB a 4.27*(R21**2.- 

Rll**2.) 

)096 



VW»4.27*(R22**2.- 

R12**2.) 

3097 



WTCMC 3 RH01*VB+RH02*VW 

)098 



WRITE(11,*)' 

OUTPUT' 

3099 



WRITEdl,*) 


)100 



WRITE(11,*)' 

BUMPER THICKNESS 3 ',T1,'CM' 

)101 



WRITEdl,*)' 

WALL THICKNESS 3 ',T2,'CM' 

3102 



WRITE(11,*)' 

MIN. WEIGHT 3 ' ,WT, ' CM ' 

J103 



WRITEdl,*)' 

CMC MIN. WEIGHT 3 ’ ,WTCMC,'KG' 

)104 



WRITE(11,*) 


3105 



WRITE(11,*) 


3106 



WRITEdl,*) 


3107 



GO TO 10 


) 108 

45 


READ(10,*)D 


)109 



READ(10,*)RH0P 


1110 



READ(10,*)S 


3111 



READ(10,*)RH0 


3112 



WRITE(11,*) ' 

MADDEN' 

3113 



WRITE(11,*) 


3114 



WRITE(11,*)' 

INPUT' 


VAX FORT® 
SAI USRDIS 



[MPACT5VM$MAIN 


17 -Aug-1987 15:27:34 

VAX FORTRAN 

'HI 




17-Aug-1987 15:19:56 

SAIJJSRDISK 

■ 

)115 


WRITE (11 ,*) 



1 

)1 16 


WR I TE (11,*) 1 

PROJECTILE DIAMETER IN CM = ' ,D 



)117 


WRITE (11 ,* ) 1 

PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 


■ 

3118 


WRITEdl,*) 1 

BUMPER/WALL DENSITY IN GM/CUBIC CM = ' ,RHO 


■ 

)119 


WRITEdl,*)’ 

BUMPER/WALL SEPARATION IN CM = ' ,S 


■1 

3120 


WRITE(11 ,*) 




3121 


T1T=0.0 



| 

3122 


T2T=0.0 



1 

J123 


DO 46 J=1 ,64 

ORIGINAL PAGE IS 


* 

J124 

1125 


K=J+8 

V=FLOAT(K) 

OF POOR QUALITY 


I 

)126 


CALL MADDEN ( V , D , RHOP , S , RHO , T1 , T2 , WT , WTCMC ) 


■ 

3127 


T1T=T1T+T1*XPV(K) 




3128 


T2T=T2T+T2*XPV(K) 



■ 

3129 

46 

CONTINUE 



1 

)130 


Tl=TlT/3 .1335 




3131 


T2=T2T/3.1335 




3132 


WT=T1+T2 



1 

3133 


R12=211. 




3134 


R22=211 .+T2 




3135 


Rll=211 .+T2+S 



i 

3136 


R21=211 .+T1+T2+S 



1 

3137 


VB=4.27*(R21**2.- 

R1 1**2 . ) 



3138 


VW=4.27*(R22**2 .- 

R 1 2**2 . ) 



3139 


WTCMC=RHO*( VB+VW) 
WRITE( 11 ,*) 1 


■ 

3140 


OUTPUT' 


m 

3141 


WRITE( 11 ,*) 




3142 


WRITE (11,*) 1 

BUMPER THICKNESS = ' ,T1 , ' CM' 


| 

3143 


WRITE (11,*) 1 

WALL THICKNESS = ',T2,'CM' 


1 

3144 


WRITEdl,*)’ 

MIN. WEIGHT = 1 ,WT, 'CM' 



3145 


WR I TE (11,*) 1 

CMC MIN. WEIGHT = ', WTCMC, 'KG' 



3146 


WRITE(11,*) 



■ 

3147 


WR I TE (11 ,*) 



w 

3148 


WRITE(11,*) 




3149 


GO TO 10 




3150 

55 

READ(10,*)D 



1 

3161 


READ(10,*)RH0P 




3152 


READ(10,*)RH01 



■i 

3153 


READ(10,*)RH02 



■ 

3154 


READ( 10 ,*)S 



1 

3155 


READ(10,*)XL2 




3156 


WRITEdl,*) 1 

WILKINSON' 


■ 

3157 


WRITE(11,*) 



■ 

3158 


WRITEdl,*) 1 

INPUT' 


JH 

3159 


WRITE(11,*) 




3160 


WRITE(11,*)' 

PROJECTILE DIAMETER IN CM = ' ,D 


I 

3161 


WRITEdl,*)' 

PROJECTILE DENSITY IN GM/CUBIC CM = ' ,RHOP 


■ 

3162 


WRITE(11,*) ' 

BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 



3163 


WRITEdl,*)' 

WALL DENSITY IN GM/CUBIC CM = ’ ,RH02 


■ 

3164 


WRITE(11,*)' 

BUMPER/WALL SEPARATION IN CM = ' ,S 
WALL MATERIAL CONSTANT = ' ,XL2 


1 

3165 


WRITE(11,*)' 


m 

3166 


WRITE(11,*) 




3167 


T1T=0.0 



l 

3168 


T2T=0.0 



1 

3169 


DO 56 J=1 ,64 




3170 


K=U+8 



i 

i 

3171 


V=FLOAT(K) 





[MPACT5 VM$MAI N 
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)172 



CALL WILKINSON( V,D,RH0P,RH01,RH02,S,XL2, 

3173 


& 

T1 ,T2 ,WT,WTCMC) 

1174 



T1T=T1T+T1*XPV(K) 

1175 



T2T=T2T+T2*XPV(K) 

1176 

b6 


CONTINUE 

1177 



T1=T1T/3.1335 

)1 78 



T2=T2T/3.1335 

)179 



WT=T1+T2 

)180 



R12=211 . 

1181 



R22=211.+T2 

1182 



Rll=211 .+T2+S 

1183 



R21*211.+T1+T2+S 

1184 



VB=4. 27* (R2 1**2. -R1 1**2. ) 

1185 



VW*4. 27* (R22**2.-R 12**2.) 

1186 



WTCMC=RH01*VB+RH02*VW 

1187 



WRITE(11,*)' OUTPUT' 

1188 



WRITE(11 ,*) 

1189 



WRITE(11,*)' BUMPER THICKNESS * '.Tl.'CM' 

WRITE ( 11 ,*) ' WALL THICKNESS = ',T2,'CM' 

1190 



1191 



WRITEdl,*)' MIN. WEIGHT = ',WT,'CM' 

1192 



WRITE(11 ,*) 1 CMC MIN. WEIGHT = '.WTCMC.'KG' 

1193 



WRITEdl,*) 

1194 



WRITE(11,*) 

1195 



WRITEdl,*) 

1196 



GO TO 10 

1197 

65 


READ(10,*)D 

1198 



READ(10,*)RH01 

1199 



READ(10,*)S 

1200 



REA0(10,*)THETA 

1201 



REAO( 10 ,*) XN 

1202 



REA0(10,*)E1 

1203 

***** 

MOOIFIED BURCH ***** 

1204 



WRITE{ 11 ,*) 1 MODIFIED BURCH' 

1205 



WRITEdl,*) 

1206 



WRITE(11,*)' INPUT' 

1207 



WRITE(11,*) 

1208 



WRITE(11,*)' PROJECTILE DIAMETER IN CM = ' ,D 

1209 



WRITE(11,*)' BUMPER DENSITY IN GM/CUBIC CM = ' ,RH01 

1210 



WRITE(11,*)' BUMPER/WALL SEPARATION IN CM = ' ,S 

1211 



WRITE(11,*)' IMPACT ANGLE FROM NORMAL IN DEGREES = '.THETA 

1212 



WRITE(11,*)' NUMBER OF PLATES TO PENETRATE AFTER FIRST', 

1213 


& 

' BUMPER = ' ,XN 

1214 



WRITEdl,*)' BUMPER YOUNGS MODULUS IN GM/CM-SQUARE' , 

121b 


& 

' SEC = ' ,E1 

1216 



WRITE{11,*) 

1217 



T1T=0.Q 

1218 



T2T=0.0 

1219 



DO 66 J»l,64 

1220 



K=J+8 

1221 



V=FLOAT(K) 

1222 



CALL BURCH (V,D,RH01,S, THETA, 

1223 


& 

XN,£1,T1,T2,WT,WTCMC) 

1224 



T1T=T1T+T1*XPV(K) 

122b 



T2T=T2T+T2*XPV(K) 

1226 

66 


CONTINUE 

1227 



T1=T1T/3.1335 

1228 



T2=T2T/3 .133b 


i 


:MPACT5VM$MAIN 
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1229 


WT=T1+T2 

1230 


R12=211 . 

1231 


R22=211 .+T2 

1232 


Rll*211.+T2+S 

)233 


R21=211 .+T1+T2+S 

)234 


VB=4 .27*( R21**2 . -R1 1**2 . ) 

1235 


VW=4.27*(R22**2.-R12**2. ) 

)236 


WTCMC=RH01*VB+2.81*VW 

1237 


WRITE(ll t *)' OUTPUT' 

1238 


WRITE(11,*) 

)239 


WRITE(11,*)' BUMPER THICKNESS = 

1240 


WRITE (11 ,*) ' WALL THICKNESS * ' 

1241 


WRITE(11,*)' MIN. WEIGHT = ' ,WT 

1242 


WRITE ( 11 ,*) 1 CMC MIN. WEIGHT * 

1243 


WRITE( 11 ,*) 

1244 


WRITEdl,*) 

1245 


WRITE (11,*) 

1246 


GO TO 10 

1247 

10 

CONTINUE 

1248 


STOP 

1249 


END 


',T1, 'CM' 
T2,'CM' 

'CM' 

,WTCMC, 'KG' 


PROGRAM SECTIONS 


Name Bytes 

0 SCODE 4826 

1 $PDATA - 743 

2 $ LOCAL 1028 

Total Space A1 located 6597 


Attributes 

PIC CON REL LCL SHR EXE RD NOWRT LONG 
PIC CON REL LCL SHR NOEXE RD NOWRT LONG 
PIC CON REL LCL NOSHR NOEXE RD WRT LONG 


iNTRY POINTS 

Address Type Name 
0-00000000 IMPACT5VMSMAIN 


/ARIABLES 


Address 

Type 

Name 

Address 

2-0000019C 

R*4 

D 

2-000001DC 

2-00000190 

1*4 

IV 


2-00000194 

1*4 

NRUNS 


** 

R*4 

R22 

2-000001E0 

2-000001B8 

R*4 

RHOP 

2-000001C4 

2-000001A8 

R*4 

T1 

** 

2-00000 1D4 

R*4 

THETA 

2-000001A4 

2-000001B0 

R*4 

WT 

2-000001B4 


Type 

Name 

Address 

Type 

Name 

R*4 

El 

2-000001A0 

R*4 

H 

1*4 

J 

** 

1*4 

K 

R*4 

Rll 

** 

R*4 

R12 

R*4 

RHO 

2-000001BC 

R*4 

RH01 

R*4 

S 

2-000001CC 

R*4 

SY1 

R*4 

TIT 

2-00000 1AC 

R*4 

T2 

R*4 

V 

** 

R*4 

VB 

R*4 

WTCMC 

2-000001C8 

R*4 

XL2 
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3001 SUBROUTINE NYSMITH(V,D,H,T1 ,T2,WT,WTCMC) I 
)002 DMAX=0 ,24*H*V**-0 . 2 ™ 
J003 IF(D.GT.DMAX)THEN 

)004 WRITE(11,*)' NO SOLUTION— PROJ . OIA. TOO LARGE FOR NYSMITH' ■ 
)005 GO TO 16 I 
3006 END IF 

J007 T1=1.93*V**0.18*D**1.91/H**0.91 m 
)008 T2 a l .86*T1 ■ 
)009 16 CONTINUE ■ 

3010 RETURN 

3011 END ■ 


>ROGRAM SECTIONS 



Name 

Bytes 

Attributes 







0 

$CODE 

150 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

1 

SPDATA 

49 

PIC CON REL 

LCL 

SHR 

NOEXE 

RD 

NOWRT 

LONG 

2 

SLOCAL 

8 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 


Total Space Al located 

207 









-1NTRY POINTS 

Address Type Name 
0-00000000 NYSMITH 


/ARIABLES 


Address Type 

Name 

Address 

Type 

Name 

Address Type 

Name 

AP -000000080 R*4 

D 


R*4 

DMAX 

AP-OOOOOOOCO R*4 

H 

AP -000000 140 R*4 

T2 

AP-000000040 

R*4 

V 

AP-000000180 R*4 

WT 


.ABELS 

Address Label 


0-00000095 


16 


D001 
3002 
)003 
3004 
)005 
)006 
)007 
)008 
)009 
)010 
)011 
)012 
)013 
3014 
)015 
JO 1 6 
501 7 
JO i 8 
)019 
J020 
J021 
J022 
J023 
J024 
3025 
J026 
)027 
)028 
)029 
)030 
)03i 
1032 
)033 
)034 
)035 
J036 
J037 
J038 
J039 
J040 
JU41 
J042 
J043 
J044 
)04b 

3046 

3047 

3048 

3049 

3050 

3051 

3052 

3053 

3054 

3055 

3056 

3057 


17-Aug-1987 15:27:34 VAX FORTRAfi 

17-Aug-1987 15:19:56 SAIJJSRDISK 

SUBROUTINE B0EING(V,D,RH0P,RH01,RH02,S,XL2,SY1 ,SY2, THETA, 

& XN,E1 ,T1 ,T2,WT,WTCMC) 

***** PEN4 ***** 

Tl=0.16 

V=V*3280. 

0=0/30.48 

RP=D/2.0 

RH0P=RH0P*1 .94 ORIGINAL PAGE IS 

RH01=RH01*1 .94 OF POOR QUALITY 

RH02=RH02*1 .94 ' J 

NITSP=0 

NITSP=NITSP+1 

NP1=0 

TlP=Tl/30.48 

T2P=FT2P ( RHOP , V , RP , SY1 , THETA , RH01 ,SY2,D, RH02 .TIP) 

WT=T1P+T2P 

IF(NITSP.EQ.1)THEN 

T1P1=1.1*T1P 

T2P1=FT2P(RH0P, V.RP.S Yl, THETA, RH01 ,SY2,D,RH02,T1P1) 

WT1=T1P1+T2P1 

ENDIF 

IF(WT1.GT.WT)THEN 

T1P1=0.82*T1P1 

T2P1=FT2P(RH0P,V,RP,SY1,THETA,RH01 ,SY2,D,RH02,T1P1) 

WT1=T1P1+T2P1 

590 IF(WT1.GT.WT)THEN 

GO TO 601 
ELSE 

T1P=T1P1 

T2P=T2P1 

WT=WT1 

NP1=NP1+1 

IF(NP1.EQ.100)THEN 

WRITE(11,*) 1 NO CONVERGENCE IN PEN4‘ 

GO TO 557 
ENDIF 

T1P1=0.9*T1P1 

T2P1=FT2P(RH0P,V ,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 

WT1=T1P1+T2P1 
GO TO 590 
ENDIF 
ELSE 

579 T1P=T1P1 

T2P=T2P1 
WT=WT1 

T1P1=1.1*T1P1 

T2P1=FT2P ( RHOP, V,RP,SY1, THETA, RH01,SY2,D,RH02, TIPI) 

WT1=T1P1+T2P1 
IF(WT1.GT.WT)THEN 
GO TO 601 
ELSE 

NP1=NP1+1 

IF(NP1.EQ. 100) THEN 

WRITE(11 ,*) 1 NO CONVERGENCE IN PEN4‘ 

GO TO 557 
ENDIF 
GO TO 579 



■OEING 


1058 

3059 

)06U 

1061 

)062 

)0b3 

)U64 

)065 

3U66 

)067 

3068 

3069 

3070 
J071 

3072 

3073 

3074 

3075 
)U76 
)077 

3078 

3079 

3080 

3081 

3082 

3083 

3084 

3085 

3086 

3087 

3088 

3089 

3090 

3091 
)092 
J093 
3094 
)095 
3096 
30 9 7 
3098 
30 99 

3100 

3101 

3102 

3103 

3104 

3105 

3106 

3107 

3108 

3109 

3110 

3111 

3112 

3113 

3114 
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END IF 
ENDIF 

601 CONTINUE 
D a 30. 48*0 
RHQP-RHOP/1.94 
RH01 a RH01/l .94 
RH02 a RH02/l .94 
TlP a 30.48*TlP 
T2P a 30.48*T2P 
IF{TlP/D.LE.0.4)VF a 4100 
IF(TlP/D.GT.0.4)VF a 4986*(TlP/D)**0.21 
VF a VF+4000. 

IF(V.LE.VF)TH£N 

WRITE(11,*)' INSIDE OF PEN4 LIMITS' 

Tl a TlP 
T2 a T2P 
GO TO 1102 
END IF 

557 CONTINUE 
***** WILKINSON ***** 

V a V/3280. 

Tl-O. 604*0**2. *RH0P/S 
Tl a Tl*SQRT(V/(XL2*RH01*RH02) ) 

RATIO a D*RHOP/ (T1*RH01) 

IF ( RATIO.GT .1 .0 ) T2 a Tl 
IF(RATI0.LE.I .0)T2 a Tl/RATI0 
***** . MODIFIED BURCH ***** 

VB*V*3280. 

DB-D/2.54 

CM a SQRT(El/RH01) 

CM-CM/30.48 

SB a S/2.54 

IF(TH£TA.LE .0.001 )G0 TO 125 
CHI a TAN( THETA) -0.5 

F2 a 0.5-1 .87*(T1B/D)+(5.*T1B/D-1 .6)*CHI**3.0 
F2 a F2+(l .7-12.*TlB/D)*CHI 
F3 a 0.32*(TlB/D)**0.83 

F3 a F3+0.48*(TlB/D)**0.33*(SIN(THETA))**3.0 
T2F a D*( (F1+0.63*F2)/XN)*(CM/V)**2.29 
T2F a T2F*( D/S) **0.71 
T2N a F3*(CM/V)**l .33*D/XN 
IF(T2N.GE.T2F)T2B a T2N 
IF(T2N.LT.T2F)T2B a T2F 
T2B a T2B*2.54 
IF(T2B.GT.T2)NREGI0N=3 
IF(T2B.GT.T2)T2 a T2B 
GO TO 155 
125 CONTINUE 

TlBl a Tl/2.54 

NITSB=0 

XKl a (DB/XN)**1.71*(CM/VB)**2 .29/SB**0 .71 

VDELTA a 0.0 

DELTA3 a 0.52 

1099 0ELTA2*2 .33* ( 1 .-1 .57*DELTA3 ) 

DELTAl a l .33*(2.*DELTA3-1 . ) 

VOELTAl- ( 1 ,/DELTAl ) **DELTA1* ( 2 .8*XK1/ ( DELTA2*DB**0 . 57 ) ) **DELTA2 , 
VDELTAl a VDELTAl*(l .58*XK1*DB**0.57/DELTA3)**DELTA3 


BOEING 


3115 

3116 

3117 

3118 

3119 

3120 

3121 

3122 

3123 

3124 

3125 

3126 

3127 

3128 

3129 

3130 

3131 

3132 

3133 

3134 

3135 

3136 

3137 

3138 

3139 

3140 

3141 

3142 

3143 
1144 
3145 

1146 

1147 

3148 

3149 

3150 

3151 

3152 

3153 

3154 
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IF(VDELTA1.LT.VDELTA)THEN 
DELTA1=1 ,33*(2.*DELTA3-1 .04) 

T1B=DELTA1*VDELTA 
T2B=VDELTA-T1B 
GO TO 499 
ENOIF 

VDELTA=VDELTA1 
DELTA3 a DELTA3+0 .02 
IF(DELTA3.GT.0.63)THEN 
T1B=DELTA1*VDELTA 
T2B=VDELTA-T1B 
GO TO 499 
ENOIF 

GO TO 1099 
499 CONTINUE 

***** COMPARISON OF MODIFIED BURCH AND WILKINSON ***** 

199 CONTINUE 

T10W=Tl/2.54 

F10W=1 .58*(DB/T10W)**0.57+2.80*(T10W/DB)**0 .57 
T2BT10W=(F10W/XN)**1.71*(CM/VB)**2.29*DB**1.71 
T2BT10W=T2BT10W/SB**0 .71 
T2BT10W=T2BT10W*2 .54 
RATI0B=(DB*RH0P)/(RH01*T1B) 

T2WT10B=0 .364*D**3 .*RHOP*V/ ( XL2*RH02*S**2 . ) 
IF(RATIOB.GT.1.0)T2WT10B=T2WT10B*RATIOB 
IF(T2BT10W.GT.T2)T2=T2BT10W 
T2B=T2B*2.54 

IF(T2WT10B.GT.T2B)T2B=T2WT10B 

T1B=T1B*2.54 

IF(T1B+T2B.LT.T1+T2)THEN 

T1*T1B 

T2=T2B 

ENOIF 

155 CONTINUE 

1102 IF(T2.LE.0.01)THEN 
SIGMA=SY2/144. 

T2=3099. 1/SIGMA 
ENOIF 

156 RETURN 
END 


VAX FORTRAN 
SAI USRDISK 


1 7-Aug- 1987 15:27:34 VAX FORT* 
17-Aug-1987 15:19:56 SAI_USRDIS 

JOOl FUNCTION FT2B(DB,T1B,XN,CM,VB,SB) 8 
)002 F1=2.42*(DB/T1B)**0 . 33+4.26* (T1B/DB) **0.33 ■ 
)003 Fl-Fl-4.18 

)004 FT2B»(F1/XN)**1.71*(CM/VB)**2.29*DB**1.71/SB**0.71 ■ 
)005 RETURN ■ 
)006 END 


PROGRAM SECTIONS 


Name 

Bytes 

Attributes 







0 $CODE 

163 

PIC CON REL 

LCL 

SHR 

EXE 

RD 

NOWRT 

LONG 

2 5L0CAL 

4 

PIC CON REL 

LCL 

NOSHR 

NOEXE 

RD 

WRT 

LONG 

Total Space Allocated 

167 









•NTRY POINTS 

Address Type Name 
0-00000000 R*4 FT2B 


/ARIABLES 

Address Type Name 

AP -000000100 R*4 CM 
AP -000000080 R*4 TIB 


Address Type Name 

AP-000000040 R*4 DB 
AP-000000140 R*4 VB 


Address Type Name 

** R*4 FI 

AP-OOOOOOOC0 R*4 XN 


)0U1 

J002 

)003 

)004 

J005 

)006 

)007 

J008 

)009 

J010 

mu 
JO 1 2 
J013 
J014 

J015 999 

J016 


17-Aug-1987 15:27:34 
17-Aug-1987 15:19:56 


FUNCTION FT2P(RH0P,V,RP,SY1, THETA, RH01,SY2,D,RH02, TIP) 
A=1.33*RH0P*(V*RP)**2. 

B=8.0*SY1*EXP(-3.125E-04*V)/C0S(THETA) 

C=1.33*RH0P*RP**2.0 
D1=RP*RH01/C0S( THETA) 

XK1=1.67*(RH0P/(2.*SY2))**0.31 
XK1=XK1*(0.281*D*RH0P/RH02)**0.33 
XK1=XK1*C0S( THETA) 

C1P1=(A-B*T1P)/(C+D1*T1P) ORIGINAL PAGE IS 

IF(C1P1.LE.0.001)THEN ORIGINAL 1 AGE IS 


FT2P=0.0 


OF POOR QUALITY 


GO TO 999 
ENDIF 


FT2P=XK1*C1P1**0.31 

RETURN 

END 


VAX FORTRAN 
SAIJJSRDISK 


PROGRAM SECTIONS 


Name 

Bytes 

0 $C0DE 

222 

2 ILOCAL 

4 

Total Space Allocated 

226 


Attributes 

PIC CON REL LCL SHR EXE RD NOWRT LONG 

PIC CON REL LCL NOSHR NOEXE RD WRT LONG 


INTRY ‘POINTS 

Address • Type Name 


Address Type Name 

** R*4 C 

AP-00000018@ R*4 RHOl 
AP-000000100 R*4 SY1 
AP -000000080 R*4 V 

.ABELS 

Address Label 


O-OOUOOOOU 

R*4 

FT2P 



• 

/ARIABLES 

Address 

Fype 

Name 

Address 

Type 

Name 

★* 

R*4 

A 

** 

R*4 

B 

AP-000000200 

R*4 

0 

** 

R*4 

D1 

AP-000000040 

R*4 

RHOP 

AP-OOOOOOOC0 

R*4 

RP 

AP-00U000280 

R*4 

TIP 

AP -000000 140 

R*4 

THETA 


0-000000D9 999 



17-Aug-1987 15:27:34 VAX FORT 

17-Aug-1987 15:19:56 SAI USRD 


)001 SUBROUTINE MADDEN(V,D,RHOP,S,RHO,Tl ,T2,WT,WTCMC) 

3002 V«V*1.E05 

)003 T1 =0 .009*SQRT ( V ) *RH0P*D**2 .0 

)004 T1=T1/(S*RH0**1.5) 

)0Ub • T2*T1 

)006 RETURN 

J007 END 


PROGRAM SECTIONS 


Name 

Bytes 

0 $CODE 

69 

Total Space Allocated 

69 


Attributes 

PIC CON REL LCL SHR EXE RD NOWRT LONG 


iNTRY POINTS 


Address Type Name 
0-00000000 MADDEN 

/ARIABLES 


Address Type 

Name 

Address Type 

Name 

Address Type 

Name 

AP-0U000Q08@ R*4 

D 

AP-00000014@ R*4 

RHO 

AP-0000000C0 R*4 

RHOP 

AP-00000018@ R*4 
AP-UOOOOO240 R*4 

T1 

WTCMC 

AP-0000001C0 R*4 
• 

T2 

AP-000000040 R*4 

V 


UNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTH$SQRT 


ORIGINAL PAGE IS 
OF POOR QUALITY 


17-Aug-1987 15:27:34 VAX FORTRAN 

17-Aug-1987 15:19:56 SAIJJSRDISK 

1001 SUBROUTINE WILKINSON( V,D,RH0P,RH01 ,RH02,S,XL2, 

)002 & T1,T2,WT,WTCMC) 

1003 T1=0. 604*0**2. *RHOP/S 

)004 Tl a Tl*SQRT ( V/ ( XL2*RH01*RH02 ) ) 

1005 RATI0=D*RH0P/(T1*RH01) 

)006 IF(RATI0.GT.1.0)T2=T1 

1007 IF(RATIO.LE .1 .U)T2=T1/RATI0 

1008 RETURN 

1009 END 


Bytes Attributes 

87 PIC CON REL LCL SHR EXE RD NOWRT LONG 
87 

;NTRY POINTS 

Address Type Name 
0-00000000 WILKINSON 


>ROGRAM SECTIONS 
Name 
0 $CODE 

Total Space Allocated 


/ARIABLES 

Address Type Name 

AP-000000080 R*4 D 
AP-OOOOOOOCO R*4 RHOP 
AP -000000040 R*4 V 


Address Type Name 

** R*4 RATIO 

AP-000000180 R*4 S 
AP-000000280 R*4 WT 


Address Type Name 

AP-000000100 R*4 RH01 
AP-000000200 R*4 T1 
AP-0000002C0 R*4 WTCMC 


: UNCTIONS AND SUBROUTINES REFERENCED 
Type Name 
R*4 MTHSSQRT 



1 7 -Aug- 1987 15:27:34 
17-Aug-1987 15:19:56 


)001 
J002 
)UU3 
)0U4 
)005 
)006 
3007 
)008 
)009 
)U10 
>011 
J012 
)013 
JO 14 
3015 
)016 
)017 
)0i8 
)019 
)020 
)021 

3022 

3023 

3024 

3025 

3026 

3027 

3028 

3029 

3030 

3031 
1032 

3033 

3034 
)035 

3036 

3037 

3038 

3039 

3040 

3041 

3042 

3043 

3044 

3045 

3046 

3047 

3048 

3049 

3050 

3051 


***** MODIFIED BURCH ***** 

SUBROUTINE BURCH( V , D , RH01 , S , THETA , 

& XN,E1,T1,T2,WT,WTCMC) 

VB»V*3280. 

DB*D/2.54 

CM*SQRT(E1/RH01) 

CM=CM/30.48 

SB-S/2.54 

IF(THETA.LE .0.001) GO TO 425 
CHI*TAN( THETA) -0.5 

F2=0.5-1.87*(T1B/D)+(5.*T1B/D-1.6)*CHI**3.0 
F2=*F2+(1 .7-12.*TlB/D)*CHI 
F3*0 . 32*( T1 B/D ) **0 .83 

F3=F3+0. 48*(T1B/D)**0. 33* (SIN( THETA) )**3.0 
T2F a 0*( ( Fl+O .63*F2 )/XN )*( CM/ V ) **2 .29 
T2F=T2F*( D/S) **0.71 
T2N«F3*(CM/V)**1.33*D/XN 
IF(T2N.GE.T2F)T2B*T2N 
IF(T2N.LT.T2F)T2B=T2F 
T2B 3 T2B*2.54 
IF(T2B.GT.T2)NREGI0N=3 
IF(T2B.GT.T2)T2»T2B 
GO TO 499 
425 CONTINUE 

NITSB-0 

XK1»(DB/XN)**1.71*(CM/VB)**2.29/SB**0.71 

VDELTA»0.0 

0ELTA3=0.52 

1099 DELTA2*2 .33*( 1 .-1 .57*DELTA3 ) 

D£LTA1»1 ,33*(2.*DELTA3-1. ) 

VDELTA1* ( 1 ./DELTA1 ) **DELTA1*( 2 ,8*XK1/ (DELTA2*DB**0 .57 ) )**DELTA2 
VDELTA1=VDELTA1*(1.58*XK1*DB**0.57/DELTA3)**DELTA3 
IF( VDELTA1 .LT.VDELTA)THEN 
DELTA1=1.33*(2.*DELTA3-1.04) 

T1-DELTA1*VDELTA 
T2-VDELTA-T1 
GO TO 499 
ENDIF 

VDELTA-VDELTA1 
DELTA3*DELTA3+0.02 
IF(DELTA3.GT. 0.63) THEN 
T1-0ELTA1*VDELTA 
T2»V0ELTA-T1 
GO TO 499 
ENDIF 

GO TO 1099 
499 CONTINUE 

Tl a Tl*2.54 

T2=T2*2.54 

RETURN 

END 


VAX FORTRA 
SAI USRDfl 


